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Abstract: Due to its light mass of order Aqcd, the strange quark can play a spe- 
cial role in Chiral Symmetry Breaking (xSB): differences in the pattern of xSB in 
^ ! the limits Nf = 2 {muymd 0, physical) and Nf = 3 {mu,md,ms 0) may 

5^ I arise due to vacuum fluctuations of ss pairs, related to the violation of the Zweig 

rule in the scalar sector and encoded in particular in the O(p^) low-energy constants 
L4 and Lq. In case of large fluctuations, we show that the customary treatment of 
SU (3) X SU (3) chiral expansions generate instabihties upsetting their convergence. 
We develop a systematic program to cure these instabilities by resumming nonper- 
turbatively vacuum fluctuations of ss pairs, in order to extract information about 
xSB from experimental observations even in the presence of large fluctuations. We 
advocate a Bayesian framework for treating the uncertainties due to the higher or- 
ders. As an application, we present a three-flavour analysis of the low-energy tttt 
scattering and show that the recent experimental data imply a lower bound on the 
quark mass ratio 2ms/ {f^u + ^d) > 14 at 95% confidence level. We outline how 
additional information may be incorporated to further constrain the pattern of xSB 
in the Nf = 3 chiral limit. 

Keywords: [Spontaneous Symmetry Breaking, QUD, (Jhiral Lagrangians. 
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1. Introduction 



Light quarks have their own hierarchy of masses. On one hand, 771^ and are 
much smaller than any intrinsic QCD scale, and their non-zero values induce only 
small corrections to the SU{2) x SU{2) chiral limit, in which = rrid = 0. A 
systematic expansion in rriu and rrid, keeping all remaining quark masses at their 
physical values, defines the two-flavour Chiral Perturbation Theory (xPT) ^ . On the 
other hand, the mass of the strange quark is considerably higher (see e.g., Ref. 
and references therein for recent determinations); indeed, it is nearly of the order of 
Aqcd, the characteristic scale describing the running of the QCD effective coupling. 
Nevertheless, a slowly convergent SU{3) x SU{3) chiral expansion is conceivable [Q]. 
This is suggested from a comparison of the kaon mass Mk with the mass scale 
Ah ~ 1 GeV of (strange) QCD bound states not protected by chiral symmetry: 
Mj^/Aj^ ~ 0.25 — 0.30. Due to the rather speciflc value of m^, the strange quark 
plays a special role among all six quarks: 

i) rris is small enough to be used as an expansion parameter (at least in some 
restricted sense) and to relate properties of QCD vacuum in the SU (3) x SU (3) 
chiral symmetry limit = rrid = = to observable quantities. 

ii) Unlike mu,md, the strange quark mass is sufficiently large, rUg ~ Aqcd, to 
infiuence the magnitude of order parameters characteristic of the SU (2) x SU{2) 
chiral limit m„, rud = with fixed at its physical value. 

iii) At the same time, is not large enough to suppress loop effects of massive ss 
vacuum pairs. This is to be contrasted with heavy quarks Q = c,b,t for which 
rriQ ^ Aqcd and the effect of QQ pairs on the vacuum structure is expected 
to be tiny. 

The above remarks single out the role of massive strange sea quarks, and suggest 
a possibly different behaviour for Nf = 2 and Nf = 3 chiral dynamics. The origin of 
this difference clearly appears in connection with the possibility that in the vacuum 
channel {J^'-' = 0"'"+) the Zweig rule and the 1/Nc expansion break down. This is 
strongly suggested by scalar meson spectroscopy sum rule studies 0, ^, |^, |[, as 
well as by instanton- inspired models ||^^. Furthermore, the enhancement of Zweig- 
rule violating effects of ss pairs on chiral order parameters has a natural theoretical 
interpretation as a consequence of fluctuations of the lowest eigenvalues of Euclidean 
QCD Dirac operator, in particular of their density ||10[- These fluctuations would only 



affect quantities dominated by the infrared end of the Dirac spectrum |]rT|; however, 
this is precisely the case of chiral order parameters such as the quark condensate 
and the pion decay constant. (Most of observables not protected by chiral symmetry 



^It should also be visible in fully unquenched lattice simulations. 
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are not especially sensitive to small Dirac eigenvalues and they have no particular 
reason to break the Zweig rule or the 1/Nc expansion.) Fluctuations of small Dirac 
eigenvalues lead to a large long-range correlation between 0"*" massive ss and massless 
uu + dd pairs. This correlation enhances the SU (2) x SU (2) order parameters 

S(2) = - lim (uM)|m,=physical (1-1) 
F(2)^ = lim -F^|m^=physical (1-2) 

by a contribution which is induced from vacuum ss pairs. The "induced condensate" 



and "induced decay constant" ITT], [T^ are proportional to and vanish in the 
SU{?>) X SU{'i) chiral limit = = = 0. As a result the Nf = 3 condensate 
11(3) and the decay constant F(3)^ can be substantially suppressed compared to the 
corresponding two-flavour order parameters, 

S(2) > S(3) = S(2)U^=o (1.3) 
F{2f > F{3f = F(2)V.=o . (1.4) 

The existence of this paramagnetic effect and its sign can be expected on general 



theoretical grounds [Tl|, but its magnitude depends on the size of fluctuations of small 
Dirac eigenvalues, which is hard to infer from first principles. A general discussion 
of the interplay between chiral order and fluctuations in the QCD vacuum can be 
found in Ref . [|1^ . 



The main question to be asked is how can the effect of vacuum fluctuations on 
chiral symmetry breaking be detected experimentally. Recall that two-flavour order 
parameters are most easily accessible via low-energy vrvr scattering. Using accurate 
recent data Jl^, we have inferred values for the Nf = 2 condensate and decay 



constant; expressed in suitable physical units, we found ITH] 



A-l^)- '"":";!"' " =°-81±0.07 (1.5) 

n n 

Z{2)= =0.89 ±0.03. (1.6) 

TT 

The fact that both X{2) and Z{2) are rather close to one indicates that, as long as 
rus is kept at its physical value, the effect of nonzero mu-,md is indeed small. This 
in turn suggests that the standard two-flavour xPT is a well-behaved expansion [||]; 
its leading order, described by the decay constant = F(2)^ and by the Nj = 2 
quark condensate S(2) = F^B, is dominant. On the other hand, the three- flavour 
order parameters S(3) and -F(3)^ are more difficult to pin down, since they require 
an extrapolation to = 0. The latter necessitates the use of three-flavor xPT, 
including more observables such as the masses and decay constants of the whole 
octet of Goldstone bosons, the K — tt form factors, K — tt scattering amplitude, etc. 
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The Nf = 3 xPT involves more low-energy constants starting in order O(p^), and 
higher orders are likely to be more important than in the two-flavour case. Most 



existing analyses P, O, 18] are based on an explicit assumption that the effect of 



vacuum fluctuations of ss pairs on order parameters is small: it is usually assumed 
that the two parameters of the O(p^) Lagrangian Fq = F{3) and S(3) = FqBq are 
such that Fq ^ and (m^ + m(i)E(3) ^ F^M^, i.e., that the corrections due to 
nonvanishing rris can be treated as a small perturbation. 

A closely related assumption concerns the smallness of the two O(p^) LEC's 
Lg(/i) and L^fi) which describe the \a.Yge-Nc suppressed and Zweig-rule violating 
effects of fluctuations in the vacuum channel. Experimental information on the actual 
size of these two constants has been rather scarce; for long time it was customary to 
posit LliMp) = (-0.2 ±0.3) x lO'^ and Ll{Mp) = (-0.3 ±0.5) x lO'^ as an input to 
both one- loop |T6[ and two-loop calculations [1^ . More recently, attempts of 
indirect estimates of L4 and Lg have appeared, all pointing towards a small positive 
values compared to the old Zweig-rule based estimates mentioned above. Rapidly 
convergent sum rules for the correlator {{Uu){ss)) |^, |^, [^, ^ yield a rough estimate 
Lq{Mp) = (0.6 ± 0.3) X 10"^, while from the analysis of it — K sum-rules |jT9| it has 
been concluded that U^{Mp) = (0.2 ± 0.3) x 10^^. The last conclusion has been 
confirmed in a recent two-loop fit to the scalar form-factors [^. The point is that 
the effect of such small shifts on order parameters is amplified by large coefficients: 
with the above estimates, S(3) and F{3Y can be suppressed compared to E(2) and 
F(2)^ respectively by as much as a factor of 2. In this way, vacuum fluctuations of 
ss pairs could lead to a particular type of instability in three-flavour xPT. 

The main purpose of the present work is to investigate instabilities in Nf = 3 
xPT that would specifically arise from the (partial) suppression of order parameters 
S(3) and F^, and to propose a systematic nonperturbative modification (resumma- 
tion) of the standard xPT recipe that could solve the problem. We assume that 
the whole expansion of relevant observables in powers of rus is globally - though 
slowly and at most asymptotically - convergent. The problem may occur with the 
enhancement of particular terms of the type m^Lg or msL^ that appear with large co- 
efficients and can be identified as arising from fluctuation of vacuum ss pairs. These 
terms are responsible for the important "induced contributions" to S(2) and -F(2)^, 
explaining why S(2) and F(2)^ at physical could be substantially larger than 
their = limit S(3) and -F(3)^ respectively. We show that in order to solve this 
particular problem it is not necessary to modify the standard chiral counting rules 
as in generalized xPT [^. The modification we propose is more modest: within 



the standard expansion scheme in powers of quark masses and external momenta, 
it appears sufficient to resum the fluctuation terms driven by rrisL^ and ttIsLq in 
the usual perturbative reexpression of order parameters msS(3) and F^ in terms of 
observables such as M^, M^, ■ ■ ■ and physical decay constants. This resummation is 
of importance for the purpose of extracting the value of Nf = 3 order parameters 
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from experiment. 

The possible effects of vacuum fluctuations in three-flavour xPT and their re- 
summation are discussed in Sections 2 and 3. These two sections are focused on 
Goldstone boson masses and decay constants, which are the observables directly en- 
tering the reexpression of order parameters mE(3) and Fq. In our approach, the 
influence of higher xPT orders (O(p^) and higher) is encoded into a few parame- 
ters referred to as "NNLO remainders" , which are kept through the whole analysis 
whatever their values. The latter depend on the model one takes for the higher or- 
der counterterms and one hopes they remain reasonably small independently of the 
model used. The result of this part of our article is an exact expression of ^4,-^5, Lq 
and Lg in terms of three fundamental parameters 

and four NNLO remainders. Using these expressions inside the xPT formulae for 
various additional observables, one can hope to pin down the values of X(3), Z{3) and 
r for a given set of assumptions about higher orders (NNLO remainders). The logical 
structure of the problem naturally calls for a Bayesian statistical type approach [^]. 
As a first application we consider the three- flavour analysis of vrvr scattering, since 



today rather accurate data exist in this case and we know from past studies |TT], |T3 
that a strong correlation exists between the value of r and the characteristics of 
the two-flavour chiral limit as revealed in low-energy vrvr scattering. In Section 4 
a quantitative analysis of this correlation is presented for the first time. On the 
other hand, the Nf = 3 order parameters X{3) and Z{3) cannot be extracted from 
the TCTT data alone. In Section 5 we survey some possibilities of learning about 
these fundamental order parameters from n — K scattering, rj ^ in decays, OPE 
condensates and sum rules and, last but not least, from lattice simulations with three 
fully dynamical fermions: we present the corresponding extrapolation formulae using 
our resummed xPT formulation. 



2. Convergence and instabilities of Nf = 3 chiral expansion 

We first recall the general structure of three- flavour xPT P] , emphasising where and 



how the instabilities due to vacuum fluctuations of ss pairs could possibly show 
up. Unless stated otherwise, a typical quantity subject to the expansion in powers 
of running quark masses m„, m^, rUs will be thought of as a connected QCD correla- 
tion function of quark currents (V, A, S, P) with external momenta fixed somewhere 
in a low-energy region away from the singularities generated by Goldstone bosons. 
We will take as a working hypothesis that the usual low-energy observables, e.g., 
Goldstone boson masses, decay constants, form factors and scattering amplitudes 
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(at particular kinematical points), when linearly expressed through such QCD corre- 
lation functions exhibit optimal convergence properties. While a similar assumption 
is implicitly made in the standard off-shell formulation of xPT Q , we will shortly 
argue that in the presence of important vacuum fluctuations this assumption should 
be understood as a restriction: observables that are not linearly expressible in terms 
of QCD correlators, e.g., ratios of Goldstone boson masses, need not admit a well 
convergent perturbative treatment and they should be treated with a particular 
care. This selects for instance F^M^, F^ and F^F^Aj^k (where At,k denotes the 
vri^'-scattering amphtude), but rules out M^. 

2.1 The bare xPT series 

The chiral expansion of symmetry-breaking observables in terms of the three lightest 
quark masses mu,mrf, is actually not a genuine power series expansion, due to 
the presence of chiral logarithms, which reflect infrared singularities characteristic 
of the chiral limit. One can nevertheless give an unambiguous scale- independent 
meaning to the renormalized coefficients of each power of individual quark masses. 
An observable A can be represented as a formal series 

A = ^mim'}m[Ajki[mu,md,ms; Bo, Fo; L[{fi) . . . L\q{^)]C{{^i) . . .C^o(/^);- • •]> 

(2.1) 

where the coefficients Ajki are defined in terms of the constants contained in the 
effective Lagrangian: i) The basic order parameters Bq and Fq which are related to 
the three-flavour chiral limit of the quark condensate and of the pion decay constant 
respectively, 

S(3) = - hm iuu), Fo = F(3)= lim F^ , b^ = -^; (2.2) 



zz; the 10 O(p^) LEG'S L[(/i), zzzj the 90 0(/) LEC's C'^ifi) and eventually 
higher-order counterterms. All LEC's are renormalized at the scale /i. In addition, 
the Ajki depend logarithmically on the quark masses through the Goldstone boson 
masses in the loops, and this dependence is such that for each jkl the coefficient Aj^i 
is independent of the scale /i. The representation ( |2.1| ) has been explicitly worked 
out for some observables to one and two loops |T^, |TB|, ^ and there is no doubt 



that it extends to all orders of the chiral expansion. We shall refer to the expansion 
expressed in the form ( p.l|) as a bare expansion, to emphasize that it is entirely 
written in terms of the parameters of the effective Lagrangian - no reexpression of 
the latter in terms of observable quantities has been performed. It is crucial that 
even before one starts rewriting and reordering the series (|2.1| ) in powers of Goldstone 
boson masses, the full renormalisation of the bare expansion ( p.l[ ) can be performed 
order by order in quark masses. Consequently, the coefficients Aj^i are finite as well 
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as cut-off and renormalisation scale-independent for all values of quark masses and 
of (renormalized) LEC's in the effective Lagrangian. 

In view of possible applications, we are concerned with practical questions related 
to the convergence properties of the bare xPT expansion (|2.1| ) in QCD. The latter 
will depend on the values of running quark masses and on the values of the LEC's at 
the typical hadronic scale ~ Mp set by the masses of non-Goldstone hadrons. In 
particular, one should question the convergence of the bare chiral expansion for the 
actual values of quark masses and not just in the infinitesimal vicinity of the chiral 
limit. In the real world, all three quarks uds are sufficiently light, 

muiAn), mdiAn) < ^.(Ah) < Ah , (2.3) 
to expect a priori some (at least asymptotic) convergence of the three-flavour bare 



xPT series. As pointed out in Refs. |0, |TT|, instabilities of the latter can nevertheless 



occur due to fluctuations of massive ss pairs in the vacuum. The importance of such 
pairs is measured by the strength of the effective QCD coupling; i.e., comparing 
with Aqcdi rather than with the hadronic scale A^- Furthermore, the impact of 
these fluctuations is proportional to m^. Hence, instabilities due to fluctuations of 
vacuum quark-antiquark pairs turn out to be particularly relevant for strange quarks 
and could manifest themselves when two- and three-flavour chiral expansions are 
compared. 

It has been argued 0, |TT1 that fluctuations of ss pairs lead to a partial suppression 



of the three- flavour condensate S(3), reducing the relative importance of the first 
term in the bare expansion of the Goldstone boson masses. We can consider for 
instance the Ward identity related to the mass of the pion (from now on we neglect 
isospin breaking and take = = m): 

FImI = 2mS(3) + 2mmsZ' (2.4) 



327r2 M2 ' M| 



The parameters Z'^ and A are defined in terms of the LEC's L^lfi) , L^ifi) and log- 
arithms of Goldstone boson masses (their expression is recalled in App. |^). Vac- 
uum fluctuations of ss-pairs show up in the term rngZ^. For the physical value of 
rris ~ Aqcd, the corresponding O^p"^) term 2mm gZ^ can be as important f^, 0, H 
as the leading-order condensate term 2mE(3). Even then, the remainder c?^, which 
collects all O(p^) and higher contributions, can still be small: d^^ <^ 1. In other 
words, vacuum fluctuations need not affect the overall convergence of the bare chi- 
ral expansions such as ( |2.1| ) or ( p.4| ) at least for some well-defined selected class of 
observables. 



-7- 



2.2 The role of NNLO remainders 

Let us write a generic bare expansion in a concise form 



A = A 



LO 



A 



NLO 



AM. 



(2.5) 



Eq. (|2.5| ) is an identity: Alo collects leading powers in quark masses in the bare 
expansion ( p.l| ) (e.g., the condensate term in Eq. (|2.4|)), A^lo consists of all next-to- 
leading contributions (the second and third terms in Eq. ( p.4| )), whereas A5A stands 
for the sum of all remaining terms starting with the next-to-next-to-leading order 
(NNLO). In Eq. (|^, the latter is denoted as 6{F^M^) = d„. 

With this setting, A can be identified with the exact (experimental) value of the 
observable A. Usually, Alo corresponds to the 0{p^) contribution, Anlo to 0{p'^) 
and ASA collects all higher orders starting with 0{p^) ^. 6A will be referred to as 
"NNLO remainder". A precise definition of 6A involves some convention in writing 
the argument in the chiral logarithms contained in the one- loop expression for Anlo- 
To illustrate this point consider the typical next-to-leading expression: 



A 



NLO 



mqUlqi 



qq' 



O-qq' 



(/i) + ^ajyA;pQ(/i) 
PQ 



or 



A 



NLO 



E 



iria 



PQ 



(2.6) 



(2.7) 



corresponding respectively to a leading-order term Alo = O(mquark) and Alo = 
0(1). Here g, g' = {u,d,s) and P,Q label Goldstone bosons. We have introduced 
the loop factor in the general case of unequal masses M: 



1 M| log(M| //i^) - M| log(M2 



327r2 



Ml 



(2. 



which in the limit of equal masses becomes: 

1 



kppijj,) 



327r2 



log 



Ml 



+ 1 



:2.9) 



In Eqs. (p^.6D-(pr7|), the constants aqqi^jj) [aq(/i)] are expressed in terms of 0{j>'^) 
LEC's L\{ij) multiplied by the appropriate powers of Bq and Fq. These constants 
are defined in the chiral limit and are consequently independent of quark masses, 
similarly to the known numerical coefficients a^? ['^^'^l- 

The only requirement from xPT is that Anlo reproduces the 0{p^) behaviour 
in the limit of small quark masses mquark 0; i.e., when the Goldstone boson 

^The case of a quantity whose expansion only starts at 0(p*) or higher, requires particular care. 
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masses Mp in the loop factors ( |2.8|) and ( p.9|) are replaced by their respective leading- 
order contributions. Once this mathematical condition is satisfied, different ways of 
writing the arguments of the chiral logarithms for physical values of quark masses 
merely correspond to different conventions in defining the NNLO remainders SA. For 
observables of the form of Eqs. ( p.6| )-( P7D at O(p^), we will use the convention which 
consists in writing in Eq. ( |2.8| ) the physical values of the Goldstone boson masses 
Mp] alternatively, we could have used the sum of LO and NLO contributions to Mp. 
This concerns, in particular, the expansion of the Goldstone boson masses and decay 
constants. In the latter case one has P = Q and the convention simply amounts to 
writing the O(p^) tadpoles, in the notation of Ref. 0, as 

The same rule can be applied to the unitarity corrections arising in the bare expansion 



of subtraction constants that define form factors and low-energy rni |24| and nK [19 



25| , ^] amplitudes. Such a prescription (detailed in Sec. ^TT]) will suffice for the 



quantities considered in this article. 

Not much is known about the size of NNLO remainders despite the fact that com- 



plete SU (3) X SU{3) two-loop calculations do exist for some observables [|T^, [T8|, ^ 
and the general structure of the generating functional is known to this order [|23[ . Fol- 
lowing this line, the bare expansion can be pushed further and the NNLO remainder 
SA can be represented as 

ASA = A^M + ^iM + A,1ee(/^) + • • • , (2.11) 

where the ellipsis stands for 0{p^) and higher contributions. The splitting of the 
0(jP) part ||23[ into the genuine two-loop contribution A2L (containing only O(p^) 
vertices), the one- loop contribution An (with the insertion of a single 0{p'^) ver- 
tex) and the tree 0{p^) contribution Atree depends on the renormalisation scale and 
scheme. 

Several ingredients are actually needed to estimate 5 A from the representation 
( P^.lll ). The first two terms (loop contributions) depend respectively on 0{p'^) pa- 



rameters mBojirLsBo, Fq and on O(p^) LEC's L\{^). Furthermore, the tree-level 
counterterms A^gg(/i) are built up from the 90 LEC's Cl{ji) that define the Oijf') 
effective Lagrangian. Even if some of them can presumably be determined from 
the momentum dependence of form factors, decay distributions and scattering am- 
plitudes (e.g., quadratic slopes), the remaining unknown O(p^) constants, which 
merely describe the higher-order dependence on quark masses, are probably much 
more numerous than the observables that one can hope to measure experimentally. 
At this stage some models (resonance saturation, large N^., NJL . . . ) and/or lattice 
determinations are required P^, but the large number of terms contributing to a 
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given A^gg makes the resulting uncertainty in 6A delicate to estimate. Finally, it is 
worth stressing that only the sum of the three components shown in Eq. ( |2.11|) is 
meaningful. An estimate of the size of NNLO remainders is therefore not possible 
without a precise knowledge of the O(p^) and O(p^) constants tuBq, Fq and the Lj's. 

In this paper, we do not address the problem of determining NNLO remainders 
on the basis of Eq. ( |2.11| ). We are going to show that interesting nonperturbative 
conclusions can be reached, even if we do not decompose NNLO remainders and 
investigate the behaviour of the theory as a function of their size. We are primarily 
interested in the constraints imposed by experimental data on the fundamental QCD 
SU{3) X SU{3) chiral order parameters (and quark mass ratio) 

^(3)-^, (2.12) 

_ 2m5o _ X(3) 
^^^^ ~ ~ Z{3) ' ^^-^^^ 

under various theoretical assumptions on NNLO remainders (i.e., on higher xPT 
orders). A suitable approach to this problem is provided by Bayesian statistical 
inference [^. (See App. for a brief review adapted to the case of xPT.) The 
output of this analysis is presented as marginal probability distribution functions 
for the fundamental parameters (|2.12| ) and it depends not only on the experimental 
input but also on the state of our knowledge of higher xPT orders. In this approach 
the latter dependence is clearly stated and can be put under control: the analysis can 
be gradually refined if new information on the relevant NNLO remainders becomes 
available either through Eq. ( p.ll| ) or in another way. 



We start with a very simple theoretical assumption about higher orders: the bare 
chiral expansion of "good observables" as defined at the beginning of this section, is 
globally convergent. By these words, we mean that the NNLO remainder 6A in the 



identity (|2.5| ) is small compared to 1: 

M < 1 , (2.14) 

for the physical values of the quark masses and for the actual size of O(p^) and O(p^) 
parameters. On general grounds, one expects 5 A = O^m^^^^^,). In the worst case, its 
size should be 6A = O(m^) ~ (30%)^ = 0.1, but in many situations 6A turns out 
to be 0{msm) or even O(m^) and is therefore more suppressed ^. These cases are 
usually identified as a consequence of SU{2) x SU{2) low-energy theorems. (Such 
suppressions are not claimed from arguments based on the Zweig rule, since we never 
assume the latter.) However the NNLO remainders will not be neglected or used as 
small expansion parameters in the following. 

^We take as order of magnitudes 10 % for 0{m) contributions and 30 % for 0{ms) terms. 
This can be related to the typical sizes of violation for SU{2) x SU{2) and SU{2>) x SU{^) flavour 
symmetries. 
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On the other hand, no particular hierarchy will be assumed between the leading 
OijP') and next-to- leading O(p^) components of ( p.5|) . By definition, for infinitesi- 
mally small quark masses m^^, m<j one should have 

v4nlo < Alo , = ^ ~ 1 . (2.15) 



However, due to vacuum fiuctuations of qq pairs, the condition (|2.15|) can be easily 



invalidated for physical value of ~ ^qcd- as discussed in Ref. [^, the three- 
fiavour condensate S(3) in Eq. ( |2.4| ) may be of a comparable size to - or even smaller 
than - the term m^Z*, refiecting the vacuum effects of massive ss pairs. At the same 
time, vacuum fluctuations need not affect the overall convergence of the bare chiral 
expansion ( p.4|) , i.e., the condition (|2.14| ) can still hold for "good observables" such 



as F^M^. We will call conditionally convergent an observable for which 6 A <^ 1 but 
the hierarchy condition (|2.15| ) does not hold. 

2.3 Instabilities in chiral series 

Standard xPT consists of two different steps. 

1. The flrst step coincides with what has been described above as the "bare ex- 
pansion" in powers of quark masses and external momenta. The coefficients of 
this expansion are unambiguously defined in terms of parameters of the effec- 
tive Lagrangian B(),FQ,Li. . ., independently of the convergence properties of 
the bare expansion. 

2. The second step consists in rewriting the bare expansion as an expansion in 
powers of Goldstone boson masses, by eliminating order by order the quark 
masses m and and the three-fiavour order parameters S(3), F{3) in favour 
of the physical values of Goldstone boson masses Mp and decay constants Fp. 
For this aim one inverts the expansion of Goldstone boson masses: 

2mBo = Ml + 5Z ""P^P + • • • j ' ^2.16) 

where Cp contains the low-energy constants Li and the chiral logarithms. A 
similar "inverted expansion" is worked out for F^: 

F2 = F2||l + 5^4M|, + ...j , (2.17) 



and for the quark mass ratio 



ms + m M|- 
2m ~ m!" 



I + Y,c''pMl + ...]. (2.18) 
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As a result of these two steps, observables other than M^, M|-, (already used 
in Eqs. (|2.16|) , ( p.l7| ) and ( p.l8|) ) are expressed as expansions in powers of M|, and 



log Mp with their coefficients depending on the constants Lj, Cj, etc. 

We now argue that large vacuum fluctuations of ss pairs could represent a se- 
rious impediment to the second step, i.e., to the perturbative reexpression of order 
parameters. This may happen if the bare expansion ( p^.5D of Goldstone boson masses 
and decay constants is only conditionally convergent: the leading and next-to-leading 
contributions are then of comparable size Alo ~ ^nlo, despite a good global con- 
vergence 6A -C 1. Let us concentrate on the three-flavour GOR ratio X{3), defined 
in Eq. ( p.l2|) , which measures the condensate S(3) in the physical units F^M^. In 



the definition of X{3) we can replace F^M^ by its bare expansion ( p.4|) and inves- 
tigate the behaviour of X{3) in limits of small quark masses. First of all, in the 
SU{2) X SU{2) chiral limit one obtains 

hm X(3) = ^ (m, fixed), (2.19) 

m^O ^(2j 

from the definition of the two-fiavour condensate S(2) = limm^o F^ / {2m) . On 
the other hand, if both m and tend to zero, we obtain by definition: 

lim X(3) = 1 . (2.20) 

m,ms—*0 

Consequently, as long as the three-fiavour condensate S(3) is suppressed with respect 
to the two-fiavour condensate S(2) (S(3) ~ S(2)/2 is suggested by the sum rule 
analysis in Refs. |, §]), the physical value of X{3) cannot be simultaneously 
close to its limiting values in both SU{2) x SU{2), Eq. (|2l9| ), and SU{3) x SU{3), 
Eq. O). 

We expect the limit "m 0, physical" to be a good approximation of 
the physical situation of very light u and d quarks, and thus the limiting value 
expressed by Eq. ( p.l9| ) to be quite close to the physical X{3). On the other hand, 
the variation of X[3) between rris = and physical may be substantial, due to 
important fiuctuations of the lowest modes of the (Euclidean) Dirac operator |I0| 



which correspond to a significant Zweig rule-violating correlation between massless 
non-strange and massive strange vacuum pairs The latter contribute to the 

two-fiavour condensate by the amount mgZ'^. (We see from Eq. ( p.4|) that S(2) = 
S(3) +msZ^ + ....) As long as S(3) is comparable to (or smaller than) the "induced 
condensate" rrisZ'^, the hierarchy condition between LO and NLO Eq. (|2.15|) will 
be violated for the bare expansion Eq. ( |2.4| ), in spite of a good global convergence 
dj^ -C 1. The inverted expansion of X{3), in which F^M^ is replaced by its expansion 
Eq. (0): 

^(2) = = l + 5^c^M| + ... , (2.21) 

1 H h . . . P 

S(3) 
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would then be invalidated by large coefficients, even for quite modest (and realistic) 
values of rus at 1 GeV around 150 MeV [0. In such a case, it is not a good idea to 
replace in higher orders of the bare expansion X{3) by 1, 2mBQ by M^, Fq by F^, 
etc. 

A comment is in order before we describe in detail the nonperturbative alterna- 
tive to eliminating order by order condensate parameters and quark masses in bare 
xPT expansions. The previous example of X(3) can be stated as a failure of the 
bare expansion of Let us remark that this is perfectly compatible with our 

assumption that the bare expansion of the QCD two-point function of axial-current 
divergences (i.e., F^M^) converges globally. Indeed, consider a generic observable A 
with its bare expansion Eq. ( |2.15|) . The latter unambiguously defines the coefficients 
of the bare expansion of 1/A: 

B = Bi^o + B^i^o + B SB , B = j, (2.22) 



in terms of those of A: 



BlO — ; -Bnlo — 7T~^ ; (2.23) 



1 o _ ^NLO 
, i^NLO T2 

I-XaY SA ^ _A 



One observes that the global convergence of A does not necessarily imply the con- 
vergence of 5 = 1/A . If the expansion of A is only conditionally convergent, i.e., 
if the relative leading-order contribution Xa is not close to 1, then 6(1/ A) need not 
be small - even in the extreme case M = 0. This explains the origin of instabilities 
and large coefficients in the inverted expansion such as ( p.l6|) , ( p.lTj ) or ( p.21| ). At 



the same time, it motivates the restriction to a linear space of "good observables" 
for which 6A <^ 1. The latter is assumed to be represented by connected QCD 
correlators, and a priori excludes nonlinear functions of them such as ratios. 



3. Constraints from Goldstone boson masses and decay con- 
stants 

The conditional convergence of FpMp and/or of Fp does not by itself bar experi- 
mental determination of three-flavour order parameters X{3) and Z{3). It just may 
prevent the use of perturbation theory in relating them to observable quantities. In 
this section a systematic nonperturbative alternative is considered in detail. 

The starting point is the standard bare expansion of FpMp and Fp for P = 
71, K, rj. (See Eq. ( p.4| ) for the pion case.) As discussed above, these particular 
combinations of masses and decay constants are expected to converge well, since 
they are linearly related to two-point functions of axial/vector currents and of their 
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divergences taken at vanishing momentum transfer. As long as Z{3) = Fq/F^ ~ 1, 
the convergence of Mp would be as good as for FpMp. If however Fp is only 
conditionally convergent (i.e., Z{3) significantly smaller than 1), the expansion of Mp 
could become unstable, in contrast with that of FpMp] the perturbative expansion 
of 1/Fp would then exhibit very poor convergence. We have no prejudice in this 
respect: the size of Z{3) as well as that of X{3) remain an open problem until they 
are inferred from the data. 

The expansion of FpMp and Fp can be written in the generic form (|2.5|) , denoting 
the corresponding NNLO remainders hj dp and ep respectively. The LO of FpMp is 
given by the condensate S(3) and the NLO contribution is fully determined by the 
standard O(p^) LEC's LQ^fi), Ls{fJ.) (and L7 in the case of the rj meson). Similarly, at 
LO Fp coincides with the order parameter Fq = F^(3) and the NLO contribution is 



given in terms of L^[ii) and L^^^ii). All necessary formulae can be found in Ref. [|I3 
Here, we follow the notation of the above reference and for the reader's convenience 
the bare expansions of FpMp and Fp are reproduced in App. ^ 

3.1 Pions and kaons 

For P = 7r,K the mass and decay-constant identities (Ward identities) consist of 
four equations that involve X{3),Z{3),r = m^/m, Le, Ls, -^4, -^5 and four NNLO 
remainders d^, dx, Ctt, ck- These identities - given in App. |^ - are exact as long as 
the remainders dp,ep are maintained in the formulae and no expansion is performed. 



As explained in Refs. [|TI| |T3[, we can combine the mass and decay constant 
identities (recalled in App. |^ to obtain two relations between the order parameters 
X, Y, Z and the fluctuation parameters p and A: 

X{3) = l-e{r) -[Y{3)fp/A-d, Z{3) = 1 - r]{r) - Y{3)X/4 - e. (3.1) 

where the functions of the quark mass ratio are^: 



and the following linear combinations of NNLO remainders arise: 

r + 1 / 2 \ 
d = -d^, - e(r) H dx , (3.3) 



r — 1 \ r — 1^ 

r + 1 / 2 \ 

e = -e^ - vir) + ck • (3.4) 

r — 1 \ r — 1/ 

The LEC's Lq and L4 enter the discussion through the combinations: 

M^ 

X = 32-^(r + 2)AL4 , p = 64-^(r + 2)ALe , (3.5) 



^In this paper, we take the following values for the Goldstone boson masses and decay constants: 
= 139.6 MeV, Mk = 493.7 MeV, Af^ = 547 MeV, = 92.4 MeV and Fk/F^ = 1.22. 
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where the scale-independent differences ALj = LKfi) — Lf^^{fi) involve the critical 
values of the LEC's defined as: 



2567r2 

1 

512^ 



log 



Ml 
1 

2r 



1 



log 



1287r2 (r- l)(r + 2) 

Mix 



lo£ 



M 



K 



Ml 



log 



M2j 
2 
9 



(3.6) 



Ml 2 , M 2 
^ + - log ^ 



3 log — ^ + log ' 



Ml 



Ml 



5127r2 (r + 2)(r-r 
The critical values of L4 and Lg are only mildly dependent on r; for r = 25, 

-0.26 • 10"^ Lf *(Mp) = -0.51 • 10"^ [r = 25]. 



(3.7) 



(3.8) 



The remaining two equations of the n — K system can be reexpressed as a relation 
between e(r) and Lg on one hand and between 77 (r) and L5 on the other hand: 



6(r) = 16--f[r(3)pAL8-d'. 

TT 

r/(r) = 8-^r(3)AL5-e'. 



(3.9) 
(3.10) 



These relations involve the combinations of the NNLO remainders d' = d — d^^ and 
e' = e — e-jr- At large values of r (> 15), e(r) and ?7(r) are suppressed and d' <ti d d^, 
e' <^ e ~ e^. 

The LEC's arise in Eqs. ( |3.9| ) and ( |3.10| ) through the differences: 



5127r2 
1 



, Ml 2, 

log ^ + - log ^ 



5127r2(r - 1) 
1 



M2 

31og^ + log^ 



Ml 



Ml 



M' 



K 



(3.11) 



2567r2 
1 



log^ + 21og^ 



2567r2(r 



M2 



m; 



.ft' 



Ml 
Ml 



(3.12) 



These differences combine the (renormalized and quark-mass independent) constants 
Lg, L5 and chiral logarithms so that they are independent of the renormalisation scale 
II. For r = 25, we obtain 



AL5 = Ll{Mp) + 0.67 ■ 10"^ , ALg = Ll{Mp) + 0.20 ■ 10" 



[r = 25] . (3.13) 
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3.2 Perturbative reexpression of order parameters 

The four exact equations Eqs. (lOP and (p.9D-( p?T0| ) can be used to illustrate explicitly 
the instabilities which may arise in the perturbative expression of X(3) and Z{3) in 
powers of Mp. In the perturbative treatment of three-flavour xPT [||, one uses the 
fact that Y{3) = 1 + 0{Mp) to set systematically Y{3) = 1 whenever it appears 
in the NLO term. One first uses Eqs ( p.9|) and (|3.10| ) to eliminate and 
r = uis/m in terms of AL5 and ALg. The result reads: 

^.l,sM_M^,,,,.. (3.14) 

./l^2§(l.'8^il^lAI.-2AL.]....). (3.15) 

In these formulae the quark mass ratio r appearing in the expressions for ALj has 
to be replaced by its leading order value (obtained from Eq. ( ^.15| )): 

Ml 

ro = 2^-l~24. (3.16) 



Strictly speaking, Eqs. (|3.14| ) and (|3.15| ) do not get any direct contribution from the 



vacuum fluctuation of ss pairs which violate the Zweig rule and are tracked by Lq 
and L4. The situation is quite different in the case of the identities ( p.l|) for X{3) and 
Z{3), where the terms describing fluctuations are potentially dangerous. Expressing 
them perturbatively one gets: 

X(3)^l-16|2AL«-16?Ml_±MAi, + ... (3.1T) 

TT n 

Ml ^ 2Ml + Ml ^ , , 

Z{3) = 1 - 8--iAL5 - 8 ^ " AL4 + . . . (3.18) 

TT TT 

The large coefficients characteristic of the perturbative treatment of l/(-F^M^) - 
and to some extent also of 1/F^ - now become visible. Eqs. ( |3.17| ) and ( p.l8| ) lead 
numerically to: 

X(3) = 1 - 37AL8 - 950AL6 + . . . (3.19) 
Z(3) = 1 - I8AL5 - 475AL4 + . . . (3.20) 

In Eqs. (p.lTp and (|3.18| ), the main NLO contribution comes from the M|^- 
enhanced term proportional to the O(p^) Zweig-rule violating LEC's Lg and L4. If 
the latter stay close to their critical values (corresponding to Eq. ( p.8|) for r = 25), the 
NLO contributions remain small. On the other hand, even though the values of L4 
and Lq are unknown yet, dispersive relations have been used to constrain their values: 
Lq{Mp) = (0.6±0.2) ■ 10~^ based on sum rules for the correlator {{uu + dd)ss) ^, |], 
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and Ll{Mp) = (0. 2 ±0. 3) ■ 10^ from vri^ scattering data [19]. Such estimates - rather 
different from the critical values - suggest a significant violation of the Zweig rule in 
the scalar sector, and an important role for the vacuum fluctuations of ss pairs in 
the patterns of chiral symmetry breaking. 

As an exercise, for illustrative purposes, we will now use the central values of 
the above sum- rule estimates to study the convergence of Eqs. ( |3.17| )-( prTB| ). We 



actually aim in this paper at providing a framework to determine more accurately 
the size of vacuum fluctuations directly from experimental observables. If we take 
Ll{Mp) = 1.4 ■ 10-3 and Ll{Mp) = 0.9 ■ 10-^ [|, the numerical evaluation of 
Eqs. ( p.l7|) -( ^T8D leads to the decomposition: 



(3.21) 



Qty = LO + [fluct + other] + NNLO 
^(3) - = 1 - [0-82 + 0.04] + 0(/) 

TV TT 

Z{3) = = 1 - [0.34 + 0.04] + Oip') 

(^ + 1)^= 1 - [0.00 + 0.06] +0{p'). 

For each quantity, the right-hand side is the sum of the leading term (1), the NLO 
term (the sum in brackets) and higher-order terms. The NLO term is decomposed 
into its two contributions: the first one comes from the fluctuation term (proportional 
to AL4 or ALq) and the second one collects all other NLO contributions. Fluctuation 
parameters have a dramatic effect on the convergence - they are the only terms 
enhanced by a factor of M|- in Eqs. (|3.17|) and ( p.l8|) . 

It should be stressed that the instability of the perturbative expansion of X{3) 
and Z{3) does not originate from higher order terms in the expansion ( p.4|) of F^M^. 
The latter actually factorize, and can be cleanly separated from the effect of vacuum 
fluctuations. This can most easily be established in the SU{2) x SU{2) limit ( |2.19| ). 



Using a bar to indicate that a quantity is evaluated for m = and fixed 7^ 0, one 
has 

S(2) = lim ^ ^ = S(3) + 32B^msALe + S(2)rf^ , (3.22) 
m^o 2m 



where 



Consequently, 



Ale^L^(.)-^(logf + ^logf). (3.23) 

^(3) = 11 = ^^^(l-d.). (3.24) 

^0 
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We expect the effect of nonzero m to be tiny; in particular, ALg — ALg, d.„ d-,, = 
d = 0{m1) and X{3) ~ -^(3). Similarly, using the expansion of decay constants 
displayed in Appendix A, one gets 







where e^r = e = O(m^) and: 



m; 



2 



Eqs. ( p.24| ) and ( p.25| ) are exact identities, in which the whole effect of higher orders 



is gathered into the O(m^) NNLO remainders d^^ and 67^. Even if the latter are small 
(~ 10%), the expansion of X{3) and Z{3) in powers of may break down, provided 
the magnitudes of ALq and AL4 are as mentioned above. The SU{2) x SU{2) limit of 
Eqs. ( |3.17| ) and ( p.l8| ) just exhibits the first term of such an expansion. Let us stress 
that we chose the SU (2) x SU (2) limit for simplicity here, but that the factorisation 
of higher order corrections is a general result (holding even for m 7^ 0) |]13| . 



One more remark is in order, concerning the special case in which both AL4 and 



ALg are large, but satisfy: 



AL4 = 2AL6 . (3.27) 
With this particular relation between the low-energy constants, one has: 

^ ^(3) , 2mBo 1-4 . ^ 

F 3 = -J-1 = hm — / = ^ 1 . 3.28 

This would describe a situation in which the (large) vacuum fluctuations suppress 

both the condensate S(3) and the decay constant F^, i.e., partially restore the chiral 
symmetry, and yet the ratio I1{3)/Fq = Bq remains nonzero. 



3.3 Nonperturbative elimination of 0(p'') LEC's 

We have just seen that the perturbative treatment of chiral series fails if vacuum 
fluctuations of qq pairs are large, resulting in instabilities in the chiral expansions. In 
this case, the nonlinearities in Eq. (|3.1| ), relating order and fluctuation parameters, 
are crucial, and we must not linearize these relations (hence the inadequacy of a 
perturbative treatment). 

We should therefore treat Eq. ( p.l| ) without performing any approximation. Fol- 
lowing Ref. |jl3|, we can exploit Eqs. ( p.l|) to express the chiral order parameters 
X{3) and Z{3) as functions of the fluctuation parameters p and A. The ratio of 
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order parameters Y{3) is 



5. 



r(3) = , 2|l-e(r)-d| ^ 

1 - ri{r) - e + a/[1 - r/(r) - e]^ + [p - A][l - e(r) - d] 



The nonhnear character of Eqs. ( |3.1| ) resuhs in the (nonperturbative) square root. 
We see that the behaviour of Y{3) is controlled by the fluctuation parameter p — A, 
i.e., 2Lq — L4 as can be seen from Eqs. ( p.5|) . 

The perturbative treatment sketched in the previous section corresponds to lin- 
earising Eq. ( p.29|) , assuming that the fluctuation parameter p — A <^ 1. This is 
invalid if large fluctuations occur: p and/or A are then numerically of order 1, al- 
though they count as 0{p^) in the chiral limit. Eq. (|3.29|) leads to the suppression of 
F(3), which would contribute to a stabilisation of Eq. (|3?T|) by reducing the contri- 
bution proportional to the fluctuation parameters p and A. As discussed extensively 
in Ref . ||T^ , different behaviours of the fluctuation parameters can result in a rather 
varied range of patterns of chiral symmetry breaking. 

We would like to extract information about Nj = 3 chiral symmetry breaking 
from physical observables, even in the event that the perturbative expansion breaks 
down. We could proceed in the same way as in Ref. [jl3l and express as many 



quantities as possible in terms of L4 and Lq, in order to stress the role played by 
vacuum fluctuations. In the present paper, we flnd it more convenient to take as 
independent quantities the (more fundamental) chiral order parameters X{3) and 
Z{3). We should emphasize that this corresponds to a different choice from that 
adopted in the perturbative treatment of chiral series: in standard xPT, X{3), Z{3), 
r, Fk/Ft, are iteratively expressed in terms of L4, L^, Lq, L^. In contrast, we start 
by the same four identities and express nonperturbatively L^jL^jLq, Lg in terms of 
X(3), Z{3), r, Fk/Ftj] this is a sensible treatment provided that both LO and NLO 
terms are considered. 

Keeping in mind that LO and NLO contributions can have a similar size, we treat 
as exact identities the expansions of good observables in powers of quark masses, and 
exploit the mass and decay constant identities to reexpress O(p^) and O(p^) LEC's 
in terms of r, X{3), Z{3), observables quantities and NNLO remainders. This leads 
to: 

F^ = FlZ{3) , (3.30) 
2m5o = MlY{3) , (3.31) 
2m,Eo = rMlY{3) , (3.32) 

and to: 

[YimW = I^T^^^tl - X(3) - e(r) - d\ 



^The quadratic equation for Y(3) admits two solutions, but only one of them corresponds to the 
physical case (see Ref. for more detail). 
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5127r2 V 9 ^ fi^ 



[F(3)]V 1 /,,_Mi,,_M^2 



(r-l)(r + 2)5127r2 V ^ M\ 



31og-^ + log-^ , (3.33) 



^(3)^4(/^) = ^r-l^.^.\^ - ^(3) - n{r) -e] + log ^ 



r(3) /2r + l 4r + l, M 



log77f + 27-^logxf ,(3.35) 



1287r2(r + 2) V2r - 2 M|. 4r - 4 M2 



F2 yf'?^ / /lf2 M^' 

^P)^'') - l-'M + (log ^ + 2 log 



These equations are derived from Eqs. (|3.1| ), ( p.l2| ), ( p.ll| ), and they have a much 
simpler expression in terms of ALj, introduced in Eqs. ( p^) - (p.7|) and ( p.ll| )- (|3?T^ ), 
rather than Li {i = 4, 5, 6, 8): 

y'mL, = - e{r) - X(3) - d\ , (3.37) 

Y'{3)ALs = ^^[e{r) + d'], (3.38) 

ymL, = ^(^^[1 - Vir) - Zi3) - e] , (3.39) 

y(3)AL5 = ^^[r^(r) + e']. (3.40) 

n 

The above identities are useful as long as the NNLO remainders are small. The 
presence of powers of Y{3), i.e., Bq, follows from the normalisation of the scalar and 
pseudoscalar sources in Ref. these powers arise only for 0{p^) LEC's related to 
chiral symmetry breaking (two powers for Lq,Li,Ls, one for L4 and L5), and are 
absent for LEC's associated with purely derivative terms. 

Plugging these identities into xPT expansions corresponds therefore to a resum- 
mation of vacuum fluctuations, as opposed to the usual (iterative and perturbative) 
treatment of the same chiral series. We can then reexpress observables in terms of the 
three parameters of interest X{3), Z{3), r and NNLO remainders. Before describing 
how to exploit experimental information to constrain these parameters, we should 
first comment on the case of the //-meson. 
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3.4 The r]-mass and the Gell-Mann— Okubo formula 

It remains for us to discuss the mass and decay constant of t] as constrained by 
Ward identities for two-point functions of the eighth component of the axial current 
and of its divergence. This results into two additional relations (given in App. 0) 
that involve one new NLO constant Lj and two extra NNLO remainders c?^ and 
Cr^. These two identities will be used to reexpress the LEG L7 in terms of order 
parameters and quark mass ratio, and to eliminate the decay-constant F^, which is 
not directly accessible experimentally. This new discussion is closely related to the 
old question whether the remarkable accuracy of the Gell-Mann-Okubo (GO) 
formula for Goldstone bosons finds a natural explanation within xPT and what it 
says about the size of the three-fiavour condensate. 
The combination: 



D 



GO 



(3.41) 



does not receive any 0{p'^) contribution from the genuine condensate S(3). The 
ry-mass identity (|A.10|) leads to the following simple formula for Dqq, expressed in 
units FlMl: 



A 



D. 



GO 



GO 



F^Ml 



GO ■ 



(3.42) 



where ALg was defined in Eq. (|3.11| ). Similarly, the identity for ( [A.llp can be 
put into the form: 



4/^2 
F2 ^ 3 V F2 



1 + 



1 m: 



487r2 Ff 



M2 

(2r + l)log^ 



log 



Mi 



M2 



+ eGo- (3.43) 



Eqs. ( |3.42[ ) and ( |3.43| ) are exact as long as the NNLO remainders 



d, 



GO 



F^M^ ^2 /M'2 

3^^d„ - A^S^dr 



(^GO 



F2 4 Ft e 
j|e,.i|e.-| (3.44, 



F^M^"'"' ^ F^M^ 

IT n IT n 

are included. If one follows Sec. ^]2| and treats the exact formulae ( |3.42| ) and ( p.43| ) 
perturbatively, one reproduces the O(p^) expressions given in Ref. |0, as expected. 

Remarkably, the identities ( |3.42| ) and ( ^.43| ) are simpler and more transparent 
than their perturbative version, and we find them useful to make a few numerical 
estimates which may be relevant for a discussion of the GO formula. For this purpose 
we shall use the value r = tq ~ 24 of the quark mass ratio and neglect for a moment 
the NNLO remainders dco and ecoy as well as error bars related to the experimental 
inputs on masses and decay constants. For this exercise, we also disregard isospin 
breaking and electromagnetic corrections. First, the dependence of F^ on Y{3) = 
2mBo/Ml is negligibly small: 



F2 

V 
i?2 

TT 



1.651 + 0.036 -^(3) 



(3.45) 



- 21 - 



In the estimate of Ago = Dgo/F^M^, we use F^/F^ = 1.687 and find: 

Ago = 77.70 - 74.46 + 1 = 4.24 . (3.46) 

We have spht the resuh into the three contributions corresponding respectively to 1], 
K and vr, in order to emphasize the accuracy of the formula. If we drop the decay 
constants in Ago, we obtain: 

3M2 - AMI + Ml 
Ago = Tjf = 46.06 - 50.03 + 1 = -2.97 . (3.47) 

Hence, apart from a change of sign, this more familiar definition of the GO discrep- 
ancy is of a comparable magnitude as Ago- For the reasons already stressed, the 
interpretation in terms of QCD correlation functions is more straightforward when 
F|M|, is used. 

If the origin of the GO formula were to be naturally explained by the dominance 
of the O(p^) condensate term in the expansion of F^Mj,, the order of magnitude 
of the estimate ( |3.46| ) should be reproduced by Eq. ( |3.42| ) for a typical order of 
magnitude of the 0{p^) LEC's Lg and L7 without any fine tuning of their values. 
Using Eq. (|3.9|) and neglecting the NNLO remainder d', one gets: 



16-i(r - lf[Yi3)]'ALs = (r - l)2[e(r) + d'] = 22.5 [r = ro] . (3.48) 

Hence, the typical 0{p^) contribution ALg to Ago happens to be nearly one order 
of magnitude bigger than the estimate (|3.46 ): the latter can only be reproduced by 
tuning very finely the LEG Lj: 

[Y{3)f{ALs + 2L7) ~ 1.3 X 10"^ , (3.49) 

to be compared with the above estimate [y(3)]^AL8 — 1.2 x 10^^. All this of course 
does not reveal any contradiction, but it invalidates the customary "explanation" of 
the GO formula and the standard argument against a possible suppression of the 
three-flavour condensate S(3). Therefore, the fact that the GO formula is satisfied 
so well remains unexplained independently of the size of S(3) and of the vacuum 
fluctuations. The last point can be explicitly verified: the genuine condensate con- 
tribution E(3) as well as the induced condensate m^Z^ ^ which represents an 0{j)^^ 
contribution to FpMj,, both drop out of the GO combination ( |3.41| ). 



We now return to our framework: we do not assume a particular hierarchy 
between LO and NLO contributions to chiral series, and we do not neglect any 
longer the NNLO remainders (in the case of the ?7-meson, dco aiid cgo might be 
sizeable and should be kept all the way through). It is possible to use the previous 
formulae to reexpress L7 in a similar way to Eqs. (|3.33|) - (|3.36| ): 



[r(3)]% ^ 



32(r- 1)2M2 



3F2M2 - AFlMl + F^^Ml 
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— (r — 1) [e(r) + a J 



This expression should be used to reexpress nonperturbatively L7 in terms of chiral 
order parameters {F^ is given by Eq. ( |3.43| )). We can already notice that in Eq. ( |3.5U| ), 
the first contribution, which corresponds to Aco? is 5 to 10 times suppressed with 
respect to the second term (r — l)^[e(r) + d']. 



4. Three-flavour analysis of tttt scattering 

The quantities B^m, r, Fq, L4,...,Ls appear in the bare chiral series up to NLO. 
The procedure outlined above allows us to express these eight quantities in terms of 
the masses and decay constants of Goldstone bosons. Apart from the (presumably 
small) six NNLO remainders dp and ep, this leaves three unknown parameters. We 
choose these three parameters to be the order parameters X{3) and Z{3), and the 
quark mass ratio r. The remaining terms of the O(p^) chiral Lagrangian, involving 
the LEC's Li, L2, L3, Lg, and Lio, do not affect the symmetry breaking sector of 
the underlying theory in which we are mainly interested here. 

The counting of the degrees of freedom is completely analogous to the one de- 
scribed in Ref. |T^, where a global two- loop fit has been performed to the masses, 



decay constants and A'e4 form factors. In this reference, the LEC's Li, L2 and L3 were 
also included and constrained by the three experimental results corresponding to the 
two KeA form factors at threshold and to their slope. The three parameters left unde- 
termined in the fits of Ref. are the ratio of quark masses r and the O(p^) LEC's 
L4 and Lq. The unknown 0{p^) LEC's, estimated in the above reference through 
a resonance saturation assumption, introduce some theoretical uncertainty. In our 
approach, this uncertainty is included in the NNLO remainders. We keep the latter 
throughout our calculation; however, at some point of the numerical analysis we will 
have to make an educated guess as to their sizes. The main differences between our 
approach and that of Ref. [|T7| lie in our use of a nonperturbative resummation of 



instabilities, compared to the canonical two-loop perturbative elimination of O(p^) 
LEC's; in our choice of X(3),Z(3),r as the three undetermined parameters, rather 



than L4, Lq, r; and in our treatment of higher-order remainders: in Ref. they are 
computed up to 0{p^), and the additional LEC's arising at this order are estimated 
through resonance saturation. 

In order to constrain our three independent parameters, more information is 
needed. In the present section we will examine the impact of our knowledge of tttt 



scattering observables. We have previously analysed the tttt data of Ref. |Tj] in terms 
of two- flavour order parameters, allowing a rather precise determination of them 
X{2) = 0.81 ± 0.07 and Z{2) = 0.89 ± 0.03. However, the tttt scattering parameters 
are more sensitive to the two-fiavour order parameters than to the three-flavour 



- 23 - 



ones [0, 0, m]. Expanding X(2) in SU{3) x SU{3) xPT one can obtain |TT1: 



X(2)(l - 4) = ^ [1 - <0 - d - r(3)Vi] + ^^(3) , (4.1) 
r + 2 r + 2 

where /i is a (small) combination of chiral logarithms, whose precise definition is 

recalled in Eq. ( 4.24|) below; in Ref. [11| the estimate /i ~ 0.05 was obtained. 



If e(r) is not close to 1, i.e., if r larger than 15, the term in square brackets in 
Eq. ( [4.1|) is dominant. Then X{2) has only a very weak (r-suppressed) sensitivity to 
X(3). On the other hand, its value is strongly correlated with r. We expect therefore 
vrvr scattering to provide us with valuable information about the quark mass ratio r, 
but not about the Nf = 3 order parameters X{3) and Z{3). This section is devoted 
to designing a framework testing this expectation in a quantitative way. 

4.1 Low-energy tttt amplitude 

Considerable progress has been achieved recently in the understanding of tttt scat- 
tering. The solutions of the Roy equations |Q allow one to express the amplitude, 
in the whole energy domain below 800 MeV, in terms of only two parameters (e.g., 
the scalar scattering lengths, or the parameters a^T and P^-k defined below), with 
very small uncertainty. It is therefore possible to determine experimentally these 
two parameters in a model-independent way. Furthermore, at low energy the tttt 
amplitude is strongly constrained by chiral symmetry, crossing and unitarity. It can 
be expressed, up to and including terms of order {p/AnY, as: 

A^^{s\t, u) = P{s\t, u) + J{s\t, u) + 0[{p/Ah)^] (4.2) 

where P{s\t,u) is a polynomial conveniently written (following the conventions and 
notation of Ref. p^) in the form: 

Pis\t,u) 





1 f-^TTTT ( 




3 


3 




TV 


-2Mlf + 


A2 




TT 


-2Mlf + 


A4 
ire 

TT 





2Miy 

{t-2MlY+{u-2Mlf\ , (4.3) 

in terms of six subthreshold parameters a^^r, /J^r^, Ai . . . A4. J{s\t, u) collects the uni- 
tarity cuts arising from elastic tttt intermediate states. At low energy, the contribu- 
tions of KK and 7777 intermediate states are not neglected but expanded and absorbed 
into the polynomial P{s\t,u). The general form of J{s\t,u) is dictated by successive 
iterations of the unitarity condition and it is entirely determined by the first four 
subthreshold parameters up to Oijf'): 

J{s\t,u) = U{s)Jt,t,{s) 

+ [{s~u)Vit) + W{t)]J^^{t) 

+ [is - t)Viu) + Wiu)]J^M + • • ■ , (4.4) 
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where 

U,V,W are polynomials given in terms of the parameters (y-Kn, P-Kn, ^i, ^2 and the 
ellipsis stands for (known) 0{p^) contributions. (For a more explicit form, see 
Eqs. (3.18) and (3.47) of Ref. 0].) 

As a first step, the general representation ( |4.2| ) of the low-energy tttt amplitude 
can be matched with experimental phase shifts and with the solution of Roy 



equations [30| in order to determine the subthreshold parameters aTr^,/^^,^, .... This 



has been done in Ref. |jT5|, independently of any xPT expansion or predictions. 



leading to the following values: 

= 1.381 ±0.242, /?^^ = 1.081 ± 0.023 , (4.6) 

with the correlation coefficient between the two parameters = —0.14. 

At the second stage, xPT can be used to expand the subthreshold parameters 
«7r7r, /^TTTT, Ai, . . . A4 in powers of quark masses = = m and/or m^, thereby 
constraining the possible values of chiral order parameters. Notice that the expansion 
of the subthreshold parameters is expected to converge better than that of scattering 
lengths, the latter being more sensitive to small variations of the pion mass. Similarly, 
as already discussed in Sees. || and ||, one should bear in mind the possibility of 
a strong dependence of on m^: it seems therefore preferable to consider the 
expansion of F^A^^t^, i.e., the on-shell four-point function of the axial- vector current 
(rather than the scattering amplitude itself) in powers of quark masses and external 
pion momenta. In Ref. the corresponding expansion of subthreshold parameters 
a-jr-K and /^t^tt in powers of m (with fixed at its physical value) was converted into a 
determination of the two-fiavour order parameters X{2) and Z{2). Now we consider 
the "bare expansion" of F'^Mla.^T, and of the slope parameter F^/3,r7r in powers of 
m and in order to investigate directly how tightly the available tttt experimental 
information constrains the three-fiavour condensate X(3) and decay constant Fq (or 
Z(3)), and the quark mass ratio r = ms/m. 

In order to establish the "bare" SU{3) x SU{3) expansion of F^M^a^rTr and of 
F^Pt^t, we proceed as follows. We start by rewriting LO and NLO xPT contributions 
in a form similar to Eq. (|4.2|): 



F^A^^is\t,u) = F^P'-{s\t,u) + F^.r{s\t,u) + . . . . (4.7) 

P^{s\t,u) collects all LO and NLO tree and tadpole contributions and is of the 
form ( |4.3| ) with the two cubic terms omitted. The second term collects the one-loop 
contribution of order O(p^). Both are renormalized and separately depend on the 
renormalisation scale fi. The loop part reads: 



FU'(s\t,u) = s 
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+ -[it - 2M^yjUt) + iu- 2M^yjUu)] 



+ is 



u)t 
t)u 



Here Jpp and Mpp are the standard loop functions for the Goldstone boson P (see, 
e.g., ref. [0). They are related to the functions Jpp(s) through: 



Jppis) 



M 



pp\ 



t/pp(s) — 2kpp , 

s - AMI y / ^ 1 , 

— — Jpp[s) - -kpp + 

12s o 



2887r 



2 ■ 



(4.9) 
(4.10) 



At low energy and for P = K,ri, the loop functions are replaced by their expansion 
at small s: 

-2A;pp, M},p = -hpp. (4.11) 
o 



jr 

Jpp 



In these equations, we have kpp = [log(M|,//i^) + l]/327r^. 

Multiplying Eq. ( f4.4| ) by and dropping all terms beyond Oi^p"^), one should 
recover the formula ( [4.8|) . This fact can be used to work out the bare expansion of 
the subthreshold parameters contained in F^P{s\t,u). Comparing ( [4. 41) and 
leads to: 



F:P(s\Lu) = FtP^sM.u 



+Mt k. 



+ -kKK 



--[{t-2M'^Y + {u-2M'^Y]k. 



--[{s-u)t + {s-t)u] ( k^ 



1, 1 



(4.12) 



which holds for all powers of s, t, u provided that one retains just the LO and NLO 
powers in quark masses in the expression of corresponding coefficients. In agree- 
ment with the convention explained in Sec. |2.2| , we keep the physical Goldstone 
boson masses in the arguments of chiral logarithms arising from both tadpoles and 
unitarity loops. The right-hand side of Eq. ( [4. 121) is checked to be independent of 
renormalisation scale /z. The resulting "bare expansion" of the two subthreshold 
parameters which carry information on the symmetry-breaking sector reads: 



F^Mla^^ = 2mS(3) + 16m^A + 2m(8m + m,)Z" 
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Fl(3^^ = + Ami + 2(4m + m,)e 

log —5- + 2 log 



47r2 



Ml 
Ml 



+ Flep, 



(4.13) 



(4.14) 



The subthreshold parameters are thus expressed in terms of A,Z^,i,i, which are 
defined in App. ^ as scale-invariant combinations of the O(p^) LEC's and chiral 
logarithms. In order to account for NNLO and higher chiral orders, we have added 
in Eqs. ( |4.13| )- p.l4| ) the direct NNLO remainders da^jm and ep^T^T^. (Their exact role in 
our analysis will be discussed shortly.) We would like to stress that the only quantities 
we really subject to a chiral expansion are the subthreshold parameters multiplied 
by appropriate powers of and M^. The scattering amplitude as a function of 
s, t, u is given by Eq. ( [4.2|) , which holds up to and including 0{p^) accuracy with all 
singularities and threshold factors correctly described using the physical Goldstone 
boson masses. 

We have just illustrated how the bare expansion of subthreshold parameters is 
obtained in practice. The last step now consists in replacing in the bare expansion 
( [4.13D -( ^mD the parameters of the Lagrangian ttiBq, ttIsBq L^^h), L^if^), Le^fi), Ls{fi) 
by the three basic QCD parameters X{3), Z{3),r = rris/m using the identities ( ^.301 )- 
( p.36|) , which yields: 
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F(3) ,r(3)rr^(r) 



327r2Fi^^^^ ^3 



r + 2 
7 



r + 2 



6 



1 + 



r + 2 
rri{r 



d + ^^^e-2Y{3)e' 



r + 2 
1 - Z{3) 



(r- l)(r + 2) 
[{da 



M2 

(r + 2) log 
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2) log 
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(4.15) 



(r- l)(r + 2) 
2e'l. 



Ml Ml 
(2r + l)log-^ + (4r + l)log-^ 



(4.16) 



These equations relate the two observable quantities a^^.,^ and /^tttt to the three in- 
dependent parameters X(3),Z(3),r (recall that 1^(3) = X{3) / Z{3)), and contain 
the dependence on the direct remainders da,-„-K, ep^n-K as well as on the indirect ones, 
stemming from the mass and decay constant identities. The order of magnitude of 
such remainders can be estimated and will be discussed later. 
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4.2 The Bayesian approach 

The determination of the three-flavour order parameters from the experimental 
knowledge of a,r7r and P-^t, is a subtle issue, due to the nonlinear character of Eqs. (|4.15|) - 
( [4 .161) and to the fact that our three independent parameters are subject to con- 



straints (e.g., they have to be positive). Moreover, the NNLO remainders are in fact 
unknown; we can only estimate their order of magnitude. One is led to consider 
them as sources of error, but the propagation of errors is not so straightforward, pre- 
cisely because of the above-mentioned non-linearities. For these reasons the standard 
method of maximum likelihood is inadequate in this case. A convenient approach is 
provided by Bayesian analysis as described in App. ^. 

We introduce the correlated probability density function -Pexp(«7r7r, /^tttt), 



VdetC f 

^TTTr ~" «exp \ ^ f '^U Cl2 



Pe^piann, Ptttt) = — ^ — exp ( --V'^CV 



y=[ 7^T'' l c= r^^^^M, (4.17) 

I \ Cl2 C22 / 

with 1/cii = 5^(l-p2^), I/C22 = 5|(l-Pa/3), l/ci2 = ( 1 " P^/j) /Pa/3, and whcre 

5q,, 5/3 and pap are the experimental results, uncertainties and correlations 
given in Eq. ( [4.6|) . This distribution summarizes the experimental result: the two 
numbers, Oexp and /^exp, given with an accuracy specified by the matrix C . Suppose 
now that we know, independently, the true values of (a^^r, /^^^r); Eq. ( |4.171 ) can then 



be interpreted as the probability of obtaining, when performing an experiment, the 
values (oiexp, /^exp), i-e., the actual observed values, given our independent knowledge 
of (a^TT, /^tttt)- Eq- ( |4.17| ) is completely symmetric under exchange of (a^rTr, /^tttt) with 

(c^exp; /3exp) • 

This last interpretation is suitable for our problem, in which theory relates the 
subthreshold parameters a^r^ and /J^jr to the three parameters X(3),Z(3),r and to 
the remainders: 

a^^ = A[r, F(3), Z{3), 6,, 62, 63, 6,] , (4.18) 
/3^^ = B[r,Yi3),Zi3),52,6,]. (4.19) 

The relevant remainders, denoted by i = 1, 5, are defined in Table |1|; 5q and 67 
are additional remainders discussed in the following section. Therefore we can say 
that the probability of obtaining the data effectively observed, for a given choice of 
r,r(3),Z(3),<5i, is 

P(data|r, Y, Z, S,) = Pexp[^(r, Y, Z, 5,), B{r, Y, Z, <5,)]. (4.20) 

This quantity is the likelihood of the observed experimental data. Indeed, it is 
not what we are interested in: data have certainly been observed, whereas the the- 
oretical parameters are unknown to us. Instead, the probability that we need is 
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Remainder 


Definition 




Si 


d 


0.1 




e 


0.1 


S3 


e' 


0.01 




{da,-Kn - d) + Ad' 


0.03 




{ef},n7T - e) + 2e' 


0.03 


Se 


[{1 + 2/r){d^ - d) + 2d/r]/{l - d) 


0.03 


57 


[(l + 2/r)(e,-e) + 2e/r]/(l-e) 


0.03 



Table 1: Definition of NNLO remainders, ai denotes the expected order of magnitude, 
and barred remainders are taken in the Nf = 2 chiral limit. 



P(r, F, Z, 5j I data). This object can be calculated as a result of a statistical infer- 
ence using Bayes' theorem, at the price of introducing some "subjective" a priori 
knowledge of the theoretical parameters, 7r(r, Y, Z, 6i). The result of the experiment 
is viewed as an update of our previous knowledge of the theoretical parameters, 
represented by the "prior" 7r(r, "K, Z, 5,), 

P(data|r, Y, Z, 5i) ■ 7r(r, Y, Z, Si) 



P{r,Y,Z, (5i|data) 



(4.21) 



/ dr dY dZ dSi P(data|r, Y, Z, 5i) ■ 7r(r, Y, Z, Si, 

The stronger the significance of data, the weaker will be the dependence of the final 
result on our choice of the prior. 

4.3 Choice of the prior 

The prior should refiect our beliefs about the theoretical parameters before our learn- 
ing of experimental results. If we have no reason to prefer one value to any other, 
then a fiat prior should be chosen. We certainly have to implement the require- 
ment of positivity for X(3) and Z{3), so the support of the function vr will have 
a lower boundary at X{3) = 0, Z{3) = 0. A similar requirement can be applied 
to the two-fiavour order parameters X{2) and Z{2). Using their expansions in the 
SU{3) X SU{3) xPT one can obtain fTTf 

^(3) + [1 - ^(3) - <r) - d - r(3)Vi] 

Z(3) - r^{r) - e - Y{3)g,] 



X(2)(l 
Z{2){1 



r + 2 
Z{3) + -^\1 



r + 2 



(4.22) 
(4.23) 



where, /i and gi are (small) combinations of chiral logarithms 
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(4.24) 
(4.25) 
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estimated in Ref. |Tl[ as /i ~ 0.05 and gi ~ 0.07, and barred quantities refer to the 
SU{2) X SU{2) chiral limit. The positivity of X{2) and Z{2) imphes then a lower 
bound for the indirect remainders 5i and S2, 

X{2) >0^6i< = 1 - - 1"(3) Vi , (4.26) 
Z{2) >0 ^ 62< = 1 - V{r) - Y{3)gi . (4.27) 

Actually one can also establish some upper bounds for the order parameters X{3) 



and Z{3) using the so-called paramagnetic inequalities Eqs. (pr3| ) - (|1.4|) stated in 
Ref. which translate into lower bounds for the remainders 6q and 67 defined in 

Table |, 

X(3) < X(2) ^S,> 6r = 1 - ^~^'x^s){l-lf^'^' ' ^^-^^^ 

<r 7(o\ ^ A ^ Amin i 1 - ^2 " r/(r) - ^(3)^1 
m<Z{2)^5,>5, =1 z{3){l-62) • ^^-^^^ 

The ratio of order parameters Y{3) is also bounded from above, as can be seen from 
Eq. (lOQl ), 

Y{3) < F'""" = 2i^44^ • (4-30) 
1 - r]{r) - 62 

We do not assume any further knowledge about the three-flavour order parameters. 
The hypothesis of a strict convergence of chiral series, in the sense that every order of 
the expansion should be smaller than the previous one, would lead to the choice of a 
prior concentrated around the values X(3) ~ 1 and Z(3) ~ 1. In our approach this is 
not required; we allow the data to indicate whether vacuum fluctuations destabilize 
the chiral series or not. Therefore, apart from the constraints listed above, we will 
choose a flat prior for the three- flavour order parameters. In fact, this statement 
is ambiguous, since we have introduced three different quantities, X{3), Y{3) and 
Z{3), related by X{3) = Y{3)Z{3). If we restrict the problem to a flat prior in two 
parameters among [X{3),Y{3), Z{3)] = {X,Y,Z), we can get three different prior 
p.d.f.'s in (Y,Z): 

(X, Z) : n{Y, Z)dYdZ = dXdZ = ZdYdZ ^ n{Y,Z) = Z , 

{X, Y) : n{Y, Z)dYdZ = dXdY = YdYdZ n{Y, Z) = Y , (4.31) 

(r, Z) : 7r(r, Z)dYdZ = dYdZ n{Y, Z) = 1 . 

A physical argument helps us to select one of these possibilities. When Z{3) vanishes, 
chiral symmetry restoration occurs and the marginal probabilities of chiral order 
parameters must vanish. This can be obtained if the prior p.d.f. vanishes in the 



^The statement is that chiral order parameters dominated by the infrared end of the spectrum 
of the EucHdean Dirac operator should decrease as the number of massless flavours increases. The 
paramagnetic inequalities also apply to X{Nf) and Z{Nf), for Nf — 2,3. 
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limit Z{3) — > 0, which is reahzed in the case of a flat prior in {X,Z). Henceforth, 
we will focus on this case and therefore take ^{Y, Z) = Z inside the range allowed 
by the positivity and paramagnetic constraints, and 7r(F, Z) = outside. 

Inspired by Eq. (|3TT|), we restrict the range of variation of r such that < e(r) < 
1, which yields the range 

n TV TV n 

No other previous knowledge is supposed for the quark mass ratio r, so the prior is 
taken to be flat in r over this allowed domain. 

Finally, we must discuss how to treat the remainders 6i, . . . ,67. We recall that 
they are defined as the NNLO contributions in the chiral expansions of a selected 
class of observables, for which a good overall convergence is expected. Therefore, the 
remainders should be small. This expectation can be implemented in the prior by 
considering the remainders as Gaussian distributions centered at zero, with widths 
corresponding to their expected order of magnitude. Since the size of the corrections 
in the chiral series is 30% for the three- flavour expansions and 10% for the two- flavour 
ones, we will use the following rule of thumb, already introduced in footnote H: the 
typical size of remainders is (30%)^ = 0.1 for genuine 0(mg) remainders like d,e, 
30% X 10% = 0.03 for 0{msm) combination of remainders (such as da — d, ep — e: see 
below), and 0.01 for e' = 0{e/r) and d' = 0{d/r). These assumptions for the widths 
of the Gaussian distributions of the remainders are collected in the last column of 
Table 0. 

As far as Table |I] is concerned, we still have to show that the two combina- 
tions da,TVTv — d.j, and e[^^T,T, — e,r, contributing to the remainders ^4 and ^5, are of 
order 0{mms), instead of O(m^). This follows from a SU{2) x SU{2) low-energy 
theorem: consider the subthreshold parameters a-„T, and /^tttt in the Nj = 2 chiral 
limit, mu-,mct 0; from their SU{2) x SU{2) chiral expansion in powers of m ~ 



Eqs. (32)-(33) of Ref. |T5l - it is straightforward to conclude that 



lim = 1, lim /^tttt = 1. (4.33) 

m— >0 m— >0 

Then we combine SU (3) x SU (3) chiral expansions of the subthreshold parameters 
and the mass and decay constant identities for the pion, by subtracting Eq. ( |A.1| ) 
from Eq. (|4.15| ) [and Eq. ( |A.6| ) from Eq. (|4.16| )], which yields the following equalities: 



F^M^ia^^ - 1) = 12m' A + 12m'Z' - M^F(3)4m(e + 2^ (4.34) 
+ T^^(3)' ( log ^ - log TTT - ^ + FlMl{d^^^^ - d^) , 



+F^{ep,Tvn - e.) . (4.35) 
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On the right hand side of Eq. (^4.34|) [Eq. (|4.35|) ], all the terms involving O(p^) 
leg's and chiral logarithms are of 0{w?) [0{m)]. If we divide Eq. ([4.34|) by F^Ml 
[Eq. ( [4.35|) by F^] and take the Nf = 2 chiral limit m ^ 0, all the terms vanish apart 
from the NNLO remainders. The latter must therefore fulfill the following condition: 

lim {do,,n7T - dn) = 0, lim (c/?,^^ - e^) = 0, (4.36) 

m— >0 m—*0 

SO that the difference between da,TTn and d-,, [c/j^tttt and is only of order O(mms). 

Having collected all the ingredients for our choice of the prior, we can now 
perform the integration over the NNLO remainders to obtain P{r, Y, Z|data). If we 
integrate the latter with respect to two parameters, we end up with the marginal 
(posterior) probabilities: 

P(r|data) = j dYdZ P{r,Y, Z\data) (4.37) 

P(r|data) = j drdZ P{r,Y,Z\data) (4.38) 

P(Z|data) = j dYdY P{r,Y,Z\data) (4.39) 

P(X|data) = j drdYdZ5{YZ - X) P(r, F, Z|data) 

drdz]-P{r,X/Z,Z\daia). (4.40) 

Zj 

The precise expression of P and the numerical evaluation of these integrals are de- 
tailed in App. |C|. 

4.4 Discussion and results 

Even if no information from tttt scattering data is included, the results of the integra- 
tions ( [4.37|) -( ^^ ) are nontrivial, because of the interplay of the various constraints 
imposed. For instance, the prior for r, when integrated over all other variables, will 
no longer be uniform. The probability density profiles of Fig. |I] are obtained by 
replacing P(r, y, Z, 5j| data) with 7r(r, y, Z, 5j) and integrating over all but one vari- 
able. They can be interpreted as a measure of the "phase space" for each parameter 
allowed by the theoretical constraints. The significance of vrvr data will be seen in 
the modification induced with respect to such "reference" density profiles. 

In the present framework, we may make more quantitative the well-known mech- 
anism by which a low value of r implies a suppression of the quark condensate . 
In Fig. ^ we show the reference probability density profiles for X, when the quark 
mass ratio r is taken fixed at different values, r = 25, 20, 15, 10. (The ratio of order 
parameters ^(3) behaves similarly.) 

Notice that for small values of r, for which X(3) is very small, the two-flavour or- 
der parameter X(2) should also be small. This is clear from inspection of Eqs. ( 4.22| ) 
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Figure 1: The reference probability density profiles for the order parameters X{3) (full) 
and Z{3) (dashed) [top], and for the quark mass ratio r [bottom]. 



and (|4.23|) , because e(r) ~ 1. Since we know from our SU{2) x SU{2) analysis [|l5l 
of the same vrvr data that X{2) is very close to 1, we can expect that such data will 
constrain r to stay away from such small values. This is what we observe in Fig. ^ 
in which we show the marginal probability density profile for r with the inclusion of 
the TTTT experimental data. 

These data induce a substantial change in the distribution as compared to the 
reference probability density profile (restricted phase space). However, such a broad 
distribution cannot be used for a real determination of r: at most a lower bound 
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Figure 2: The reference probability density profiles for the order parameter X{3) when 
the quark mass ratio r is taken fixed at four different values. 



for r can be given in probabilistic terms: r > 14 at 95% confidence level. The most 
probable value of r lies at r = 21 — 22, very close to standard expectations. Indeed, 
the motivation for the rearrangement of the chiral expansions, namely the possible 
importance of vacuum fluctuations and their potential to destabilize the chiral series, 
is not essential in the case of r: the fluctuation parameters L4 and Lg are absent 
from the perturbative reexpression of the chiral series of r, see Eq. ( p.l5| ), which may 
thus exhibit no instability even in the case of large fluctuations. It is therefore not 
surprising that similar values of r are obtained through the perturbative expansion 
Eq. ( p. 151 ) or the nonperturbative framework presented above. Such similarity is 
actually a welcome check of the approach advocated in this paper. 

We also show in Fig. ^ the dependence of the inferred probability density profile 
on the initial choice of the prior: the full line refers to the choice of priors described 



in Sec. while the dashed-dotted one (dotted one) corresponds to a distribution 
for the remainders with half widths (double widths), as compared to Table 0. As 
expected, broader Gaussians for the NNLO remainders lead to a flatter p.d.f. for r - 
the impact of experimental information is weakened when higher orders are allowed 
to be larger. If the most probable value of r depends to some extent on the choice 
of the prior, the same is not true for the 95% lower bound, which is almost the same 
for all distributions. 



-34- 



0.07 





0.04 - 
0.03 - 



0.02 



0.01 




Figure 3: The probability density profile for r inferred from the vrvr data (thick lines) as 
compared to the reference one (thin dotted line). The different curves represent different 
choices for the prior. The full line corresponds to the choice of priors indicated in the 
previous section. The dashed-dotted line is obtained by reducing the widths for the prior 
distributions of the remainders by a factor of two, the dotted line by increasing the widths 
by a factor of two. 



A comment on the so-called "Kaplan-Manohar ambiguity" is in order here. In 
Ref. |2T|, it was shown that the sum of the 0{p'^) and O(p^) chiral Lagrangian 



remained unchanged under a shift in the quark masses: 

m — i> m + Xmrris , rris + Xrn? , (4-41) 

and a corresponding redefinition of the O(p^) LEC's Lq ^^^. This seemingly prevents 
any determination of the quark mass ratio r. However, the Kaplan-Manohar am- 
biguity also induces a shift in the 0{p^) terms included in the NNLO remainders 
(neglected in a perturbative treatment of O(p^) expansions). We have assumed a 
good overall convergence of chiral series and therefore required small NNLO remain- 
ders, which corresponds to fixing the Kaplan-Manohar ambiguity. We point out that 
this hypothesis is not specific to our framework, and is a fairly general assumption 
when dealing with chiral series. 

On the other hand not much deviation is caused by tttt data with respect to 
the reference profiles for X{3), Z{3) and Y{3), as is clear from Fig. ^ This is 
an indication of the fact that vrvr data, as was expected, are not sensitive enough to 
X(3) and Z{3). Their determination would require the inclusion of other observables, 
which we discuss now. 
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Figure 4: The probability density profiles as inferred from the vrvr data (thick lines) 
compared to the reference one, representing the phase space (thin lines). Full lines refer 
to X{3), dashed lines to Z{3) [top] and full lines to the ratio of order parameters 5^(3) 
[bottom] . 

5. Other sources of constraints on X{3) and Z{3) 

We now briefly mention some possible sources of independent experimental informa- 
tion about the three-flavour chiral order parameters X[3) and Z{3), or equivalently 
about the 0{p'^) fluctuation parameters Y'^{3)Le and Y{3)L4. The combination of 
O(p^) LEC's (AL8-I-2AL6), invariant under the Kaplan-Manohar transformation, is 
of particular interest. It appears in many observables and has remained undetermined 
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so far. Y^{3)Ls is essentially given by r through Eq. Y^{3)ALs = 0.99 ■ lO^^ 

for r = 25 (where NNLO remainders are neglected). From the positivity of X{3) and 
its relation to the fluctuation parameter Lg Eq. ( |3.1| ), one gets the following upper 
bound for Lq: 

Y\3)ALq< -Il-ll - e(r) - d] , (5.1) 

^ ' ^ - 16(r + 2) M2 ^ V ; J , v ; 

which implies that Y'^{3)ALq < 0.98 ■ 10~^ for r = 25 (with an expected uncertainty 
of 10 % from NNLO remainders). Between the limit of no fluctuations and that of 
maximal fluctuations, the combination of LEC's Y'^{3)[ALs + 2ALq\ can thus vary 
in the range: 

1.0- 10"^ < r2(3)[AL8 + 2AL6] < 3.0- 10"^ [r = 25] (5.2) 

This combination can be perturbatively related to the LEG £3 |jl[] of the Nf = 2 
chiral Lagrangian; see Eq. (11.6) in Ref. 0]. In the limit of no fluctuations, £3 is 
a positive number of 0(1); the estimate £3 = 2.9 ± 2.4 was obtained in Ref. ^ 
under the assumption of the validity of the Zweig rule. If larger vacuum fluctuations 
make Y^{3)[ALs + 2ALq] increase, £3 decreases towards negative values with a larger 
magnitude. Thus, in principle, the size of vacuum fluctuations could be investigated 
through an accurate determination of £3 obtained from vrvr-scattering parameters. 
Unfortunately, experimental information at such accuracy is not available due to a 
large uncertainty in the scattering length Oq, which is not tightly constrained by 
available KeA decay experiments ^. Therefore, we have to consider other sources of 
information in order to constrain the Nf = 3 order parameters. 

5.1 Goldstone boson scattering and decays 

In order to estimate the size of vacuum fluctuations of ss pairs, processes directly 
involving strange mesons are required. Before sketching how our method could be 
extended to the relevant processes, we want to comment on a few estimates of the 
O(p^) LEC's L4 and Lq that are available in the literature. These estimates show a 
common feature: they rely on the standard (perturbative) treatment of vacuum fluc- 
tuations, assuming that the latter are small, but they lead to L4 and Lq significantly 
larger than their critical values. These estimates should therefore be considered at 
most as valuable hints of internal inconsistency of the perturbative treatment. 

Biittiker et al. |^ have analyzed the ttK — > ttK and vrvr —>■ KK amplitudes, 
thereby obtaining the amplitudes at threshold and in the subthreshold region. Com- 

'^The available 7r7r-scattering data was analysed in Ref. yielding £3 = —18.5 ± 16.7 and 
thus suggesting an important role played by vacuum fluctuations. If the experimental information 
on the 1 = 2 channel is replaced by a theoretical constraint concerning the scalar radius of the 
pion 132, ^q|, £3 becomes a small positive number compatible with the absence of fluctuations. 
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paring these results with the 0{p'^) xPT expansion of Bernard et ah they de- 
termine the LEC's Li {i = 1,4) and the combination Lg + 2Lq: 



10^[L^ + 2L^](Mp) = 3.66 ± 1.52. (5.3) 

The large uncertainty quoted is reflected in the uncertainty in the scattering lengths 
combination + 2ao^^. This could improve considerably once experimental results 
from 71 K atoms become available. However, the analysis is based on 0{p^) xPT, in 
which the implicit assumption is made that Y{3) = 1. From our bounds Eq. ( |5.2| ), 
taking Y{3) = 1, we obtain: 

0.27 <10^[Ll + 2Ll]{Mp) <2.3 [r = 25, y(3) = 1] . (5.4) 

It is clear that the values for Lg + 2Lq given in Ref. [P3| are barely compatible 



with the assumption of 1^(3) = 1. Furthermore, the authors of Ref. |^ estimate 
lO^LKMp) = 0.50 ± 0.39, which implies lO^AL^ = 1.0 ± 0.39 and suggests signifi- 



cant violation of the Zweig rule in the scalar sector [see Eq. ( |3.20|) ]. Both remarks 



call for a comparison of ttK scattering amplitudes with a chiral expansion treating 



nonperturbatively (possibly large) vacuum fluctuations p6 . 

Recently, Bijnens and Dhonte ||2^ have calculated the vr and K scalar form 
factors at two loops in iVj^ = 3 xPT; they then fit their results to the correspond- 
ing dispersive representation based on the solution of the multi-channel Omnes- 
Muskhelishvili equations 0, In order to obtain "decent fits" for the form factors 
at zero momentum transfer, they found that the LEC's L4 and Lq had to satisfy the 
constraint 

L^(Mp) < 0.6 X LliMp) + 0.4 ■ 10"^ . (5.5) 

If we take this to imply that [2Lg — Ll\{Mp) < 0.8 • 10"'^, then we can rewrite the 
constraint in terms of the fluctuation parameters p, A introduced in Sec. ^: 

p-A<1.6 [r = 25]. (5.6) 

Recalling Eq. (|3.29| ), the convergence of the perturbative expansion of Y{3) = 
2mBo/M^ requires p — A ^ 1. Viewed in this light, the constraint of Ref. |^ is not 
surprising: it is simply the observation that the authors, within their perturbative 
analysis, cannot find a good fit to a perturbative xPT expansion in the presence of 
large vacuum fluctuations. In addition, they find that requiring the values of the 
scalar form factors at zero not deviate too much from their lowest order values leads 
to the estimate 0.3 < 10^Ll{Mp) < 0.6. Note that this result is roughly in agreement 



with that of Ref. ||3j discussed above. According to Table 2 of Ref. |£0|, this value 
of L4 leads to a suppression of Z{3) down to half of Z{2), as can be checked from 
Eq. (p.ll). All these considerations suggest sizeable vacuum fluctuations in the scalar 
sector. 
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We outline now how additional observables coming from Goldstone boson scat- 
tering and decay could be incorporated naturally in the Bayesian machinery, in order 
to constrain the size of vacuum fluctuations, or equivalently X{3) and Z{3). The 
most obvious candidate is tiK scattering. The first step consists in using analyt- 
icity, unitarity, crossing symmetry in conjunction with experimental data in order 
to constrain as much as possible the low-energy vriC-scattering amplitude [Q. The 
second stage corresponds to an analysis similar to that in the case of vrvr scattering: 
establish a dispersive representation of the amplitude with subthreshold parameters 
(i.e., subtraction constants like a,r7r and Pnir), determine the value of these param- 
eters from experiment, then include these parameters as additional observables for 
the Bayesian analysis. The expected outcome should be more stringent constraints 
on X(3) and Z{3) |6| . 



The decay ^ Svr is a second process of interest |3J] . The standard treatment of 
this decay starts with the ratio yl^_,3^/Ax, where Ar,-,^^ is the decay amplitude and 
Ak = {M'^+ — M|,o)qcd- This ratio is then perturbatively reexpressed in terms of 
Goldstone boson masses and a single O(p^) LEG L3. As stressed previously, a bare 
chiral expansion of ratios of QGD correlation functions that are only conditionally 
convergent may exhibit instabilities. To cope with possibly large vacuum fluctua- 
tions, a better starting point could be the expansion of the quantity F^F^A^^s^r, 
which is linearly related to a QGD correlation function. Such an expansion will in- 
volve more LEG's and order parameters and allow the extraction of the latter from 
the data. Since the decay rj 3tt is forbidden in the isospin limit m„ = m^, we 
must start by extending the previous discussion of Goldstone boson masses and decay 
constants to include isospin breaking. A dispersive representation of the amplitude 
y4^_^3^ must then be written to define convenient observables. The dispersion rela- 
tions are here more involved than in the case of vrvr- and vri^'-scattering: they require 
studying rjir tctt and performing subsequently an analytical continuation to the 
(physical) decay channel. The observables thus defined, related to the behaviour of 
the decay amplitude at the center of the Dalitz plot (via slope parameters), can be 
exploited in our Bayesian framework in order to constrain further three-flavour chiral 
order parameters and to extract the quark mass ratios rriu/ms and nid/ms- 

5.2 Two-point functions and sum rules 

In this paper, our aim was in particular to pin down X{3), i.e., the Nj = 3 chiral 
condensate measured in physical units, which governs the behaviour of QGD corre- 
lation functions in the limit m„ = = = 0. Related though different quantities 
arise when the high-energy limit of the same correlation functions is studied through 
the operator product expansion (OPE). Local condensates appear then, and those 
with the lowest dimension are: 

E„ = -(0|m«|0), Erf = -(0| Jrf|0) , E, = -(0|ss|0), (5.7) 
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where the physical vacuum of the theory is denoted |0) with all the quarks carrying 
their physical masses: no chiral limit is taken. 

These OPE quark condensates occur in various sum rules for two-point correla- 



tors and could thus be determined in this framework |3^. First, they arise (multiplied 
by a mass term) in the high-energy tail of the correlators as dimension-4 operators. 
(For example, see for the case of pseudoscalar densities.) Next, in some sum 
rules, normal-ordered condensates of the type ( [5.7[ ) appear through chiral Ward 
identities. For instance, in the case of the divergence of the strangeness-changing 
vector current |3^, the strange-quark mass is determined via a sum rule with 



no subtraction, but another sum rule can be written with the subtraction constant 
{nis — m^j)(Ss — Sy), providing in principle some information on the OPE quark con- 
densates. Unfortunately, the high-energy tail of the (Borel transformed) two-point 
function involved in this case has a QCD expansion which behaves quite badly and 
prevents any accurate determination. Lastly, the OPE quark condensates arise when 
factorisation is invoked to reexpress higher- dimensional four-quark operators as the 
square of qq vacuum expectation values. 

We stress that the OPE quark condensates S^, S^, Eg have a different definition 
(and thus value) from the chiral condensates that we have considered here and in 
Ref. p[: 

S(2) = - lim (0|mu|0) = lim E„ = lim , (5.8) 
S(3) = - lim {0\uu\0) 

mu,md,ms—*0 

= lim Su = lim = lim . (5-9) 

In particular, S^, S^, S., exhibit an ultraviolet divergence that must be renormalized; 
therefore, their definition and their value depend on the convention applied. It is 
possible to relate them to S(3) using Nf = 3 xPT. For instance, if we take Eq. (9.1) 
in Ref. in the isospin limit, we get: 

= = XiS) + [¥(3)]'^ [16(r + 2)AL, + A{2ALs + AH2)] + d^,^,, 

TV TV TV 

1 1 
= 1 - -e(r) + 4[r(3)]2-^A/J2 - d + -d' + , (5.10) 

X [16(r + 2) ALg + 4r(2AL8 + AH2)] + d^., 

r — 9 M"^ r 

= 1 + -^e{r) + 4r[F(3)]2-^AiJ2 -d+-d' + d^,, , (5.11) 

where NNLO remainders are denoted c?s;u,d and c/x;;s, and the 0{p^) high-energy 
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31og^ + log^) . (5.12) 



counterterm H2 arises in the combination 

1 / M| , M2 

-256vr2(r-l) I' + 

The value of such high-energy counterterms cannot be fixed by low-energy data 
only, and their presence in the chiral expansions is merely a manifestation of the 
renormalisation-scheme dependence of the OPE quark condensates. 

An interesting relation, free from high-energy counterterms, exists between the 
OPE condensates: 

^^"''^ Z^'' = 1 - e(r) - H ^c?s;„,d ^d^., . (5.13) 

r — 1 r — 1 r — 1 

Two conclusions can be drawn from this relation. First, for r larger than 15, 
Eq. ( |5.13| ) shows that Xu^d is close to [1 — e(r) — d], while we see from Eq. (|4.1| ) 
that X(2)(l — J^) equals [1 — e(r) — d] up to 1/r-suppressed corrections. T,^^ should 
thus be very close to S(2), which was expected since the u,d quarks are very light 
and the physical world is near the Nf = 2 chiral limit. 

The second conclusion is that Xg — Xu^ can hardly be obtained from such a 
relation, since X^^cl and 1 — e(r) — largely cancel. Thus, very accurate knowledge of r 
and Xu4 would be needed to determine Xg this way. More generally, the possibility of 
significant vacuum fiuctuations of ss pairs makes it difficult to relate in a quantitative 
way S(3) and the OPE quark condensates T^ud-, S^. 

Such a relation between OPE and chiral quark condensates is naturally relevant 
for the description of i^-decays. In particular, some K ^ tttc weak matrix elements 
are related to v.e.v.'s of four-quark operators in the Nf = 3 chiral limit, thanks to 



sum rules for vector-axial or scalar-pseudoscalar correlators . These sum rules are 
evaluated by splitting the integral in two energy domains: the low-energy region is 
described by experimental data, while the high-energy behaviour is obtained through 
the operator product expansion, which involves a priori OPE quark condensates. 
However, since the sum rules are evaluated in the chiral limit rriu = fn^ = 0, 
these condensates actually reduce to the Nf = 3 chiral condensate S(3). 

The dispersive estimates of K tttt matrix elements could thus be affected 
at three different stages by significant vacuum fiuctuations of ss pairs leading to 
smaller values of S(3). Firstly, extrapolating the weak matrix elements from the 
Nf = 3 chiral limit to the physical values of the u, d, s-quarks could not be carried 
out on the basis of the usual treatment and values of LEC's of Nf = 3 xPT, since 
the latter assume from the start a leading contribution from S(3). Moreover, the 
very estimate of the sum rule could be modified because of the change in the high- 
energy behaviour of the correlator in the chiral limit. The third question concerns 
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dimension-6 four-quark condensates, which appear at higher orders of OPE and are 
often related to the square of a qq condensate through factorisation, on the basis 
of large-Nc arguments. The presence of large qq fluctuations might render such an 
assumption invalid. 

5.3 Implications for lattice simulations 

In principle, the lattice should represent a particularly favourable domain to study 
how QCD at low energy depends on the light-quark masses and how this dependence 
is connected to the vacuum fluctuations of qq pairs. Recent progress has been made 
in this field. Discretisations of the Dirac operator have been discovered with highly 
desirable qualities for the simulation of light quarks. In particular, Ginsparg- Wilson 
fermions |Q do not break chiral symmetry explicitly. A second (cheaper) option 



consists in twisted-mass lattice QCD |41|, where a parametrized rotation of the mass 



matrix allows one to restore chiral symmetry partially in observables through an 



averaging procedure. Another avenue is provided by staggered fermions [^, which 
allows one to study an odd number of flavours, at the cost of introducing unwanted 
flavour degeneracies. 

Vacuum fluctuations of qq pairs are typical sea-quark effects. The fermionic 
determinant plays here an essential role, since we are interested in chiral parameters 



dominated by the infrared end of the Dirac spectrum ||T0|, |TT[ . In order to study these 
effects on the lattice, it is therefore mandatory to generate data for 3 dynamical 
flavours. For this particular purpose, one cannot rely on quenched data (with no 
dynamical quark) or on data generated with only two dynamical quarks - even 
though they can be of interest for observables relatively insensitive to the fermion 
determinant, e.g., Mp. 

We will now illustrate, by considering "bare" expansions of "good" observables 
such as F'^Ml and F^, how lattice data could probe vacuum fluctuations of qq pairs 
and how chiral extrapolations should be dealt with if the latter turn out to be signif- 
icant. We consider a slightly optimistic situation where lattice data with reasonable 
accuracy can be generated for 3 light dynamical flavours. For simplicity, we choose 
to work in the limit of degenerate strange and light quark masses. The analysis could 
be done for independent variations of the quark masses, but the attendant compli- 
cations do not add anything essential to our conclusions. Moreover, it is likely that 
realistic lattice calculations are more easily performed in this simplified situation. On 
the lattice, we denote by the common mass of the three degenerate light flavours, 
and by M| and Fl the common mass and decay constant of the eight degenerate 
Goldstone bosons. We keep m and for the physical values of the quarks and 
F^, Ml ... for the physical values of Goldstone boson observables. 



The observables F| and FlM^ have thus the form 

^2 71^2 



M' 
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2-, 



3L4(/i) + L5(^)---^log^ +FleL, (5.15) 



where the remainders diiGL are of order m\. 

Since -Bq, -fo aiid the LEC's Lj are all defined in the Nf = 3 chiral limit and are 
thus independent of the quark masses, we may use Eqs. (|3.30|) -( |3736D to eliminate 



them in favour of the real-world parameters X(3), y(3), Z(3), r and the physical 
values of the masses and couplings of the Goldstone bosons, leading to: 

FlMl = hFlMlX{?>) + ^{3[1 - X(3) - rf] + (r - l)[e(r) + d']] 

1 Ml . Ml 



= FlZ{?>) + 6 ;:^i^.'{3[l - ^(3) - e] + (r - l)[r/(r) + e']} 



r + 2 

-6 



327r 



2 



/Ml Ml\i 
-2(log-^ + 21og^)]+F^.. (5.17) 

Taking the ratio of these equations gives M| implicitly as a function of 6 = m^/m. 
The remainders d, e enter the resulting expressions multiplied by a factor of b/r; we 
will ignore them as well as the 1/r-suppressed contributions from d', e'. 

Since the chiral expansion requires small values of ttil, while present day lattice 
simulations prefer on the order of m^, it is important to ascertain if there is a 
range of variation for mj;, in which our equations may be applied and still give valu- 
able results. There are two different conditions to be met. First, as increases, the 
degenerate mass M^ increases, eventually exceeding M^, M^; this in itself is of con- 
cern, since the chiral expansion becomes unreliable as Ml approaches values typical 
of relevant hadronic resonances. In addition, the terms logarithmic in M| become 
negative and eventually the procedure described above is unstable. Therefore, we 
will restrict the range of variation oiniL so that these logarithmic terms do not con- 
tribute more than 25% of the total. This stability criterion constrains the allowed 
region in the space of parameters X(3), Z(3), r, m/,; however, this region includes 
values of interest. For example, in Fig. ^ we show the allowed region in X(3), rrii/m 
for the illustrative choice ^(3) = 0.6, r = 25. For rriL/m < 20, all values of X{3) are 
possible. As can be seen in Fig. ^ with such a constraint on rriL/m, we satisfy also 
the first requirement since M^ does not exceed 1.5 ■ M^- 

For values of Z{3) between 0.4 and 0.8 and of r between 20 and 30, we find 
that the allowed region does not qualitatively change: so long as we keep mi/m less 
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Figure 5: Stability criterion for 2'(3) = 0.6, r = 25. Inside the circle in full line [dashed 
line], the term logarithmic in Mj, contributes more than 25% to F| [F|M|]. All NNLO 
remainders are set to zero. 



than 20, then all values of X(3) are permitted according to our stability criterion. 
Consider, then, the variation of F| and of F\M\ as functions of m^, for fixed ^(3) 
and r, as shown in Figs. and |^ in the illustrative case 2^(3) = 0.6, r = 25. From 
Eqs. ( ^.16| )-( ^1^ ), the general dependence on X(3) is apparent: the X(3)-dependence 
must vanish for small h = niL/m as well as for b = mL/iTL (r + 2)/3~ 10 (apart 
from a residual dependence from the logarithms). Therefore, the region of interest 
is 10 < 6 < 20. Fortunately, even in this restricted interval, there is considerable 
dependence on X{3), especially for FlMf^/h. A good knowledge of the spectrum 
in this range would allow us to discriminate at least between the most extreme 
possibilities (X(3) close to X(2), X(3) close to 0)^. We stress that not all the 
observables are equivalent for this study: for instance, M| exhibits a much weaker 
sensitivity to vacuum fluctuations, as shown in Fig. |^. This cancellation between Lg 

first step in this direction was proposed by considering the dependence on m and of the 
Goldstone boson masses to extract the values of some low-energy constants for partially quenched 
staggered fermions |^3| . However, this was achieved by relying on the chiral expansion of Mp, with a 
perturbative reexpression of the fundamental parameters r, ^(3), Z{?>) in terms of Goldstone boson 
masses. Moreover, the 0{p^) LEC's L4 and Lg related to the vacuum fluctuations were eliminated 
through a perturbative redefinition of the parameters involved in the chiral series. As discussed in 
Sees. ^ and |3| this procedure need not be correct in the case of large vacuum fluctuations. 
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Figure 6: Mr/M^ as a function of b 
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Figure 7: Ff^Mf^/ [Ft^M.,^)'^ /b as a function of b. Each line corresponds to a different value 
of -^(3), whereas r is set to 25. All NNLO remainders are set to zero. 



and L4 in the case of the masses is quite usual [|T3| and does not mean that the effect 
is absent for all observables, as seen from Figs. |^ and |^. 
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Figure 8: FI/F^ as a function of h. Each line corresponds to a different value of -^^(3), 
whereas r is set to 25. All NNLO remainders are set to zero. 



We have not included the NNLO remainders di and cl here, but it is a straight- 
forward exercise to take them into account. We can check their size by writing them 
in the following form: 



Fle-j, = miEi + msmE2 + m^Es , 



(5.18) 



where Di and Ei do contain 0{p^) LEC's and chiral logarithms. If we neglect the 
dependence of the chiral logarithms on the quark masses when they vary from m, 
to niL, the same quantities will appear in the analogous expression for FlMldi, Flei- 



dL 
dn 



[bFlMh 62 
FlMll r2 



iFli 



62 



D1 + D2 + D3 
Di + D2/r + D^/r'^ 
El + E2 + E^ 
Ei + E2/r + E3/r^. ' 



(5.19) 



We can check easily that the first factor is of order one in each case (see for instance 
Figs. 1^ and H). Therefore, we conclude that the remainders d^^ei are expected to 
be of order m|, i.e., of order 0.1 ■ 62/r2 (i.e., 10 % or less). Once the remainders are 
included, an accurate determination of Z(3) is quite difficult because F| does not 
exhibit a strong sensitivity to the latter. But X(3) is still accessible because of the 
strong variation in the curvature of -F|M| as a function of 6. 



6. Summary and conclusion 

The large- A^c suppressed low-energy constants Ll{^) and Ivg(/i) encode fluctuations 
of vacuum ss pairs. We have analysed the influence of these fluctuations on the 
convergence of Nf = 3 xPT. 

i) For the physical value of m^, we assume a global (though possibly slow) conver- 
gence of the SU{3) X SU{3) chiral expansion applied to low-energy connected 
QCD correlation functions and to observables that are linearly related to the 
latter: typical examples are Fp, FpMp, F^A^tt^tttt , F^F];.Aj,k^t^k, etc. For 
such quantities, the bare expansion (as defined in Sec. ^ expressed in terms of 
the order parameters S(3) and Fq and in powers of quark masses mu,md,ms 
is likely characterised by relatively small higher-order corrections starting at 
NNLO. 

ii) Vacuum fluctuations of ss pairs affect in particular the next-to-leading-order 
(NLO) contributions to F^M^ and F^ through terms m^Lg and msL^, which 
appear in xPT with large coefficients. They reflect Zweig-rule violation and 
large l/N^. corrections in the scalar channel. Unless Lg and L4 are very precisely 
tuned to their critical values Lf\Mp) = -0.26-10-3, Lf^Mp) = -0.51- 10~^ 
the vacuum fluctuation NLO contribution to F^M^ and/or to F^ becomes 
of comparable size or even larger than the corresponding leading-order (LO) 
contributions given by the three- flavour condensate S(3) and by the pion decay 
constant Fq. In this case, the expansions of X{3) = 2mS(3)/F^M^ and Z(3) = 
Fq/F^ in powers of break down despite the convergence of F^M^ and F^. 
As a result, X{3) and Z{3), which measure Nf = 3 order parameters in physical 
units, could be suppressed well below one, implying an instability of Nf = 3 
XPT. 

iii) The instability due to vacuum fluctuations upsets the standard perturbative 
reexpression of the lowest order parameters E(3) and Fq as well as quark masses 
in favour of an expansion in terms of physical Goldstone boson masses, F^^ 
and Fk- Instead, using the four mass and decay constant Ward identities, we 
nonperturbatively eliminate the LEC's L4, L5, Lg, (-^^7), Lg in terms of the order 
parameters X{3) and Z{3) , the quark mass ratio r = rris/m and four NNLO 
remainders d.„,dK,eT,, ck which collect all higher-order contributions starting 
at 0{p^). This procedure amounts to an exact resummation of the standard 
perturbative reexpression of X(3), Z{3) and r and it applies even if the vacuum 
fluctuations suppress X{3) and/or Z{3). 



iv) In this way, values of the basic order parameters X{3), Z{3) and the quark mass 
ratio r can be constrained by experimental data as long as NNLO remainders 



are under control. In order to gradually incorporate theoretical conjectures 
about the order parameters X(3) and Z{3) (e.g., their positivity, paramagnetic 
inequalities), and on higher chiral orders (i.e., NNLO remainders), we propose 
to use Bayesian statistical inference. In this approach, previous knowledge of 
parameters is encoded into the prior probability distribution function. At this 
step, some degree of arbitrariness is introduced, but the dependence on the 
choice of the prior is guaranteed to be weak if data are significant enough. 

v) We have applied this procedure to the three-flavour analysis of elastic tttt scat- 
tering, in order to test the ability of the recent high-precision low-energy data 



obtained by the E865 Collaboration |14| to constrain three-flavour chiral dy- 



namics. We have shown that, for the present state of experimental precision, 
the data cannot determine the two fundamental Nf = 3 chiral order parameters 
X{3) and Z{3). However, the low-energy vrvr data is sufficient for us to put a 
quantitative lower bound on the quark mass ratio r = rris/m > 14 at the 95% 
confidence level, and to determine the corresponding probability distribution 
function. 

vi) The Bayesian machinery is suitable for incrementally including further ex- 
perimental information on low-energy observables, and seems especially well- 
adapted in this context, in view of an extension of our analysis to valence 
s-quark effects. In particular, ttK scattering appears to be rather promising, 
due to recent progress obtained through the solution of the corresponding Roy- 
Steiner equations [RH]. We plan to address this issue in future publications [EO]. 



vii) The possible instability of Nf = 3 chiral expansions, and the prescription 
examined in this paper to treat it, are also relevant in other respects. One 
example is the evaluation of K tttt weak matrix elements from sum rules 
based on OPE: in particular, the extrapolation of the OPE condensate from the 
chiral limit to physical values of quark masses would not be possible on the basis 
of the usual treatment of Nf = 3 xPT formulas. A similar remark applies to the 
study of the quark mass dependence of Goldstone boson masses on the lattice 
and to determinations of the LEC's by this means. We have shown how such a 
program could be pursued even in the presence of large fluctuations, provided 
simulations with three dynamical flavours would be available with quark masses 
as small as the physical strange quark mass. It remains to be seen whether 
such simulations, with proper control of the continuum and thermodynamic 
limits, will be feasible in the near future. 

viii) In order to detect possibly large ss fluctuations, another strategy is conceivable: 
one might calculate and analyse as many quantities as possible within the 
standard SU (3) x SU (3) xPT, up to and including two loops, using the standard 
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perturbative reexpression of low-energy parameters ||T7|, pO|. The instability of 



such xPT expansions could appear as an internal inconsistency of the result of 
corresponding fits with the assumptions underlying the perturbative treatment 
of standard xPT. 
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A. Mass and decay constant identities 

We recall for convenience the chiral expansion of Goldstone boson masses as discussed 



in Refs. ITTI, M 



F^Ml = 2mS(3) + 2m(m, + 2m)Z' + Am'^A + Am^B^L + F^M^d^, 
F^M]^ = {nis + m)S(3) + (m^ + m){ms + 2m)Z^ + (m^ + mfA 

3lL + FlMlc 



(A.l) 
(A.2) 



The case of the 77-meson is discussed in Sec. p.4| and below. The connection with the 
standard LEC's of the Nf > 3 effective Lagrangian is: 



+m(m, + m)B?,L + FlMldK, 



a = iqbi\lm- 



5127r2 
1 

512^ 



log 



log 



M\ 



Ml 2, 



Ml 
Ml 



(A.3) 
(A.4) 



The remainders FpMpdp collect all higher-order terms, starting at the next-to-next- 
to- leading order O(m^), in agreement with the definition in Ref. [|13|, but different 



from that in Ref. [§. The combination of chiral logarithms L is ||T3 

1 



327r2 



\ Ml ^ 
3 log ^ + log 



Similar expressions are derived for the decay constants: 
F^ = F^{3) + 2(m, + 2m)f + 2m^ 



+ 



167r2 



Af ^ M' 
2m5o<ilog^ + 21og 



K 



Ml 



F| = F\3) + 2(m, + 2m)e + (m, + m)^ + mBoL + F^Ck 



(A.5) 



(A.6) 



(A.7) 
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The scale-invariant constants ^ and ^ are related to the LEC's L4 and L5 as follows: 

1 



2567r2 
1 

256^ 



log 



Ml 

/i2 



, Ml 

log^ + 21og 



(A.8) 
(A.9) 



The remainders F^ep collect the NNLO terms ©(m^). Eqs. (pD-fprTOl) in 

Sec. 3^ can be obtained by combining the previous identities to eliminate the O(p^) 
LEC's 

The identities for the rj can be recast in a form reminiscent of the Gell-Mann- 
Okubo formula: 



(A.IO) 

= 4(r - l)m'{{r - 1){2Z^ + A) - BIl] + ^F^M^d^ - AF'^M'^cIk + F^M^ , 
2>F^-AFl + Fl (A. 11) 



ri " K 

_ 2mBo 
~ 167r2 



:i + 2r) log 



Ml 
Ml 



log 



K 



Ml 



The //-mass identity involves the new LEG Z'^ = IGBqL-j. We have also introduced 
the NNLO remainders dr^ and of order 0{ml). 



B. and Bayesian statistical analysis 

One of the main achievements of SU (3) x SU (3) chiral perturbation theory consists in 
providing a consistent framework that includes all the constraints of chiral symmetry 
when one analyses processes involving the Goldstone bosons vr, K, t]. This allows one 
to express observables as expansions in powers of momenta and quark masses up to 
a given order. Each order involves new low-energy constants whose values cannot 
be determined from chiral constraints, but nevertheless provide very much needed 
insight into the mechanism of chiral symmetry breaking. 

We would like therefore to pin down (or at least constrain) LEG's arising at first 
order in SU{3) x SU (3) chiral expansions - for instance, the quark condensate and 
the pseudoscalar decay constant in the limit m„, rrid, mg - by using available 
experimental information. 

We must therefore combine several sources of information: 

• Experimental measurements of observables af. Pj^'^lai) 

Their values have (possibly correlated) uncertainties which are described by 
probability distributions. In the theoretical basis that underlies the determina- 
tion of these observables no use of xPT should have been made. For instance, 
solutions of Roy equations were required in addition to experimental phase 
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shifts to extract tttt scattering parameters. As explained in Sec. |2]^, a reason- 
able choice of observables can be derived from QCD correlation functions of 
currents and densities, away from kinematic singularities; e.g., one may choose 
TTTT subthreshold parameters. 

• Relations between the observables and theoretical parameters T„: = Ai{Tn) 

We have explained how mass and decay constant identities can be used to 
eliminate some 0{p'^) LEC's in chiral expansions of observables. We have 
chosen the theoretical parameters to be X{3), Z{3), r, and remainders that 
comprise NNLO and higher order corrections. 

• Theoretical constraints/assumptions about the values ofTn. Ck{Tn) 

Chiral order parameters are constrained: for instance, vacuum stability requires 
X{3) > 0, whereas Z{3) is positive by definition. We have also theoretical 
prejudices about NNLO remainders if we require an overall convergence of 
chiral expansions: the relative contribution of NNLO remainders must remain 
small. 

It is quite easy from these elements to construct the likelihood function 

£(T„) = P(data|T„,ii) = nPfP[A(TO] , (B.l) 

j 

which is the probability of obtaining the observed data once a particular set of 
theoretical parameters is given (within the theoretical framework H). However, what 
we want is not £, but rather the reverse quantity P(T„|data, H), i.e., the probability 
of having a particular set of theoretical parameters, taking into account the data. 

This problem of "statistical inference" has a long history. A possible solution is 
provided by Bayesian analysis p2|) which relies on Bayes' theorem: 



where tt is the "prior" distribution, i.e., the probability associated with the theoretical 
parameters before the experimental results have been taken in consideration: 

7r(T0 = P(T„|iJ) = l[Ck{T^) . (B.3) 

k 

The denominator on the right-hand side of Eq. ( p.2| ) is just a normalisation factor. 
The marginal probability associated with a given theoretical parameter is obtained by 
integrating the joint probability P(T„|data, H) over all other theoretical parameters. 

Let us mention that there is some arbitrariness in each of these ingredients. We 
could have chosen other observables, such as combinations of scattering lengths in 
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which one-loop chiral logarithms cancel ||3^. We could have added all O(p^) LEC's 
to the set of theoretical parameters, and kept "bare" chiral expansions as relations 
among the observables. Finally, we could have computed 0{p^) contributions to 
NNLO remainders and used resonance saturation to estimate the size of the remain- 



ders, following the procedure in Ref. [|T^ . 

In the present paper we have advocated a particular choice of observables as 
constituting a sensible starting point for a Bayesian analysis of data from tttc and 
ttK scattering. Other choices of prior p.d.f 's for the theoretical parameters (especially 
for the NNLO remainders) can be considered, as long as they are well motivated. 
However, the posterior probabilities should not be strongly sensitive to the choice of 
priors when a sufficient amount of experimental data is included in the analysis. 



C. Integration procedure for the analysis of tttt scattering 

In Sec. ^, we apply Bayesian methods to perform a three-flavour analysis of vrvr scat- 
tering, constructing a joint probability P(r, Y, Z, 5|data). This gives the probability 
of having a particular choice of quark mass ratio r, order parameters Y{3) and Z{3) 
and NNLO remainders 5i=i...7, once vrvr scattering data is taken into account. By 
integrating over NNLO remainders, we obtain the joint probability 

P(r, F, Z|data) = f Yld6i P{r, Y, Z, (5|data) (C.l) 
^ 1=1 

oc n(Y, Z) e{r - 0(r2 - r) 0[Z(3)] ^[F(3)] 

d5i I d62 / dd^ / d6-/ / dd^ d5i d5^ 



X 



J]G(5„a,)exp (-^V^CV^ ^[^---^(3)], (C.2) 



where oc means "equals, up to a (numerical) normalization coefficient", C is the 
error matrix for the experimental data, Eq. ([4.17|) , G(5j, ai) is the Gaussian function 
of 5i with width aj, and V is defined in terms of the chiral series for aT,.,, and /Stttt, 
Eqs. (TO)-(CT): 

V{r,Y^Zj)={^^'^Y^l'''-A . (C.3) 
y B{r, Y, Z, 6) - /5exp / 

Before any numerical evaluation, we can analytically compute some of the inte- 
grals. 64 and 65 have the same width = and appear only in the experimental 
distribution obtained from vrvr scattering. We can therefore diagonalize the latter 



R = RR^ = / , RC R^ 



Ci 
C2 



(C.4) 
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to perform the integration ^ over 64 and 65: 



j S, 5, G(54, ^4)^(55, cr4) exp (^-iv^CV^ oc exp 



1 + Cial ' 



>V = i?-V(r,F,Z,5i,2,3;54 = 55 = 0) 



(C.5) 
(C.6) 



The integrals over 5^ and ^7 are simply Gaussians integrated over semi-infinite ranges, 
and are thus given in terms of the error function Erf. We obtain finally: 



P(r, y, Z|data) oc 7r(F, Z) e{r - n) e{r2 - r) e[Z{i)] e[F(3)] 



/ rmin \ 

l_Erf (^) 



1 - Erf 



\V2cxj) 



9[Y' 



Y{3)] 



(C.7) 



In order to obtain the marginal probability for either r, X{3), Y{3) or Z{3), 
we must perform a numerical integration over 3 remainders to obtain the joint 
probability P(r, F, Z| data), and then integrate the result over two of the remain- 
ing three parameters. We restrict the integration over the remainders (5i,2,3 to the 
range [— 5(Ti,5crj] (the upper bound can be smaller for i = 1,2 due to the positivity 
constraints Eqs. ( [4.26| )-( P:.27] )). The range of integration for the two parameters that 
remain to be integrated is given by the theoretical constraints discussed in Sec. ^. 

In order to appreciate the impact of data in the Bayesian framework, it is quite 
interesting to consider the marginal probabilities obtained when no experimental 
information is included. In our particular case, this amounts to replacing the experi- 
mental distribution Pexp(«, P) by 1, or equivalently to setting the matrix C to 0. We 
obtain then the following p.d.f: 

Po(r, y, Z|data) oc 7r{Y, Z) e{r - ri) e{r2 - r) e[Z{3)] e[Y{3)] 

„5max ^^max 3 

d6, J d62 J dds n a,)e[Y^-'' - Y{3)] 



X 



1=1 



X 



1 - Erf 



\V2aJ 



/ rmin \ 



Po is just the normalized prior p.d.f, and corresponds to the "phase space" allowed by 
the theoretical constraints and assumptions on the various parameters. The resulting 
marginal probabilities can be found in Sec. ^, where they are compared with the ones 
that include experimental knowledge about vrvr scattering. 



similar procedure can be followed in the case of different widths for the Gaussian prior p.d.f. 's 
for 64 and ^5. 
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